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Abstract 

In this review, we present and discuss the most recent results on inclusive and exclusive diffrac- 
tive processes at HERA and Tevatron colliders. Measurements from fixed target experiments at 
HERMES and Jefferson laboratory are also reviewed. The complementarity of all these results is 
analyzed in the context of perturbative QCD and new challenging issues in nucleon tomography are 
studied. A first understanding of how partons are localized in the nucleon to build orbital momenta 
can be addressed with these experimental results. Some prospects are outlined for new measure- 
ments in fixed target kinematic, at Jefferson laboratory and CERN, at COMPASS, or at the LHC. 
Of special interest is the exclusive (coherent) production of Higgs boson and heavy objects at the 
LHC. Based on the present knowledge, some perspectives are presented in this direction. 
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1 Introduction 



Understanding the fundamental structure of matter requires an understanding of how quarks and gluons 
are assembled to form hadrons. Of course, only when partons are the relevant degrees of freedom of the 
processes, which we design in the following as perturbative processes. The arrangement of quarks and 
gluons inside nucleons can be probed by accelerating electrons, hadrons or nuclei to precisely controlled 
energies, smashing them into a target nucleus and examining the final products. Two kinds of reactions 
can be considered. 

The first one consists in low momentum transfer processes with particles that are hardly affected 
in direction or energy by the scattering process. They provide a low resolution image of the structure, 
which allows to map the static, overall properties of the proton (or neutron), such as shapes, sizes, 
and response to externally applied forces. This is the domain of form factors. They depend on the 
three-momentum transfer to the system. The Fourier transformation of form factors provides a direct 
information on the spatial distribution of charges in the nucleon. 

A second type of reaction is designed to measure the population of the constituents as a function 
of momentum, momentum distributions, through deep inelastic scattering (DIS). It comes from higher 
energy processes with particles that have scored a near-direct hit on a parton inside the nucleon, 
providing a higher resolution probe of the nucleon structure. Such hard scattering events typically 
arise via electron-quark interactions or quark- ant iquark annihilation processes. Nucleon can then be 
pictured as a large and ever-changing number of partons having appropriate distributions of momentum 
and spin. 

Many experiments in the world located at DESY (Hamburg), Jefferson Lab or JLab (Virginia), 
Brookhaven (New York), Fermilab (Batavia) and CERN (Geneva) can measure these processes. Both 
approaches described above are complementary, but bear some drawbacks. The form factor measure- 
ments do not yield any information about the underlying dynamics of the system such as the momenta 
of the constituents, whereas the momentum distributions do not give any information on the spatial 
location of the constituents. In fact, more complete information about the microscopic structure lies in 
the correlation between momenta and transverse degrees of freedom. New results in this direction are 
presented in this review and the complementarity of these measurements, from all experiments listed 
above, is discussed. 

2 Basics of diffraction at HERA and Tevatron 

HERA was a collider where electrons or positrons of 27.6 GeV collided with protons of 920 GeV, 
corresponding to a center of mass energy of about 300 GeV. One of the most important experimental 
results from the DESY collider HERA is the observation of a significant fraction, around 10%, of large 
rapidity gap events in deep inelastic scattering (DIS) [H [21 El H]- In these events, the target proton 
emerges in the final state with a loss of a very small fraction (xp) of its energy-momentum. 

In Fig. [11(a), we present this event topology, 7*p — X p', where the virtual photon 7* probes the 
proton structure and originates from the electron. Then, the final hadronic state X and the scattered 
proton are well separated in space (or rapidity) and a gap in rapidity can be observed in the event 
with no particle produced between X and the scattered proton. In the standard QCD description of 
DIS, such events are not expected in such an abundance since large gaps are exponentially suppressed 
due to color strings formed between the proton remnant and scattered partons (see Fig. Hl^b)). The 
theoretical description of such processes, also called diffractive processes, is challenging since it must 
combine perturbative QCD effects of hard scattering with non perturbative phenomena of rapidity gap 
formation. The name diffraction in high-energy particle physics originates from the analogy between 
optics and nuclear high-energy scattering. In the Born approximation the equation for hadron-hadron 
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Figure 1: Parton model diagrams for deep inelastic diffractive (a) and inclusive (b) scattering observed 
at lepton-proton collider HERA. The variable f3 is the momentum fraction of the struck quark with 
respect to P — P', and the Bjorken variable xbj its momentum fraction with respect to P. 



elastic scattering amplitude can be derived from the scattering of a plane wave passing through and 
around an absorbing disk, resulting in an optic-like diffraction pattern for hadron scattering. The 
quantum numbers of the initial beam particles are conserved during the reaction and then the diffractive 
system is well separated in rapidity from the scattered hadron. 

The early discovery of large rapidity gap events at HERA [1] has led to a renaissance of the physics 
of diffractive scattering in an entirely new domain, in which the large momentum transfer provides a 
hard scale. This observation has then revived the rapidity gap physics with hard triggers, as large-p_L 
jets, at the proton-antiproton collider Tevatron (see Fig. E]). The Tevatron is a pp collider located 
close to Chicago at Fermilab, USA. It is presently the collider with the highest center-of-mass energy 
of about 2 TeV. Two main experiments are located around the ring, D0 and CDF. 
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Figure 2: Schematic diagrams of topologies representative of hard diffractive processes studied by the 
proton-antiproton collider Tevatron. 



In the single diffractive dissociation process in proton-proton scattering, pp — > Xp, at least one 
of the beam hadrons emerges intact from the collision, having lost only a small fraction of its energy 
and gained only a small transverse momentum. In the analogous process involving virtual photons, 
7*p — s> Xp, an exchanged photon of virtuality dissociates through its interaction with the proton at 
a squared four momentum transfer t to produce a hadronic system X with mass Mx- The fractional 
longitudinal momentum loss of the proton during the interaction is denoted xip, while the fraction of 
this momentum carried by the struck quark is denoted /?. These variables are related to Bjorken x by 
X = P X]p (see Fig. [3]). 

Using the standard vocable, the vacuum/colorless exchange involved in the diffractive interaction is 
called Pomeron in this review. Whether the existence of such hard scales makes the diffractive processes 
tractable within perturbative QCD or not has been a subject of intense theoretical and experimental 
research during the past decade. 
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Figure 3: Diffractive kinematics. 

3 Observation of diffractive events at HERA 

3.1 The rapidity gap events 

Let us start by giving a real example of a diffractive event in HERA experiments. See Fig. HJ which 
is the (exact) experimental reproduction of Fig. [H A typical DIS event as shown in the upper plot of 
Fig. m is ep — >■ eX where electron and jets are produced in the final state. The electron is scattered in 
the backward detectoi0 (right of the figure) whereas some hadronic activity is present in the forward 
region of the detector. The proton is thus completely destroyed and the interaction leads to jets and 
proton remnants directly observable in the detector. 

The fact that much energy is observed in the forward region is due to color exchange between the 
scattered jet and the proton remnants. However, for events that we have called diffractive, the situation 
is completely different. Such events appear like the one shown in the bottom of Fig. HI The electron is 
still present in the backward detector, there is still some hadronic activity (jets) in the LAr calorimeter, 
but no energy above noise level is deposited in the forward part of the detectors. In other words, there 
is no color exchange between the proton and the produced jets. The reaction can then be written 
as ep — >■ epX. This is also called a Large Rapidity Gap (LRG) event, and constitutes an efficient 
experimental method to tag diffractive events. 

3.2 Proton tagging 

A second experimental technique to detect diffractive events is to tag the outgoing proton. The idea is 
then to detect directly the intact proton in the final state. The proton loses a small fraction of its energy 
and is thus scattered at very small angle with respect to the beam direction. Some special detectors 
called roman pots can be used to detect the protons close to the beam. 

The basic idea is simple. The roman pot detectors are located far away from the interaction point 
and can move close to the beam, when the beam is stable, to detect protons scattered at vary small 
angles. 

The inconvenience is that the kinematical reach of those detectors is much smaller than with the 
rapidity gap method. On the other hand, the advantage is that it gives a clear signal of diffraction 
since it measures the diffracted proton directly. A scheme of a roman pot detector as it is used by the 
HI or ZEUS experiment is shown in Fig. 13.21 The beam is the horizontal line at the upper part of 
the figure. The detector is located in the pot itself and can move closer to the beam when the beam is 
stable enough (during the injection period, the detectors are protected in the home position). 

^At HERA, the backward (resp. forward) directions are defined as the direction of the outgoing electron (resp. proton). 
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Figure 4: Usual (top) and diffractive (bottom) events in the HI experiment at HERA. For a diffractive 
event, no hadronic activity is visible in the proton fragmentation region, as the proton remains intact 
in the diffractive process. On the contrary, for a standard DIS event, the proton is destroyed in the 
reaction and the flow of hadronic clusters is clearly visible in the proton fragmentation region (+z 
direction, i.e. forward part of the detector). 



The third method used at HERA mainly by the ZEUS experiment is based on the fact that there is 
a different behavior in logM^, where Mx is the total invariant mass produced in the event, either for 
diffractive or non diffractive events. For diffractive events dadijf/dMx = {s / Mx)°'~^ = const, if a ~ 1 
(which is the case for diffractive events). The ZEUS collaboration performs some fits of the da/dMl 
distribution: 



as illustrated in Fig. 13.31 The usual non diffractive events are exponentially suppressed at high values of 
Mx- The difference between the observed da/dM^ data and the exponential suppressed distribution is 
the diffractive contribution. 



3.3 The Mx method 



da 



D + cexp(61ogM|.) 
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Figure 5: Scheme of a roman pot detector. The roman pot detectors are located far away from the 
interaction point and can move close to the beam, when the beam is stable, to detect protons scattered 
at vary small angles. 
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Figure 6: Illustration of the Mx method used by the ZEUS collaboration to define difFractive events. 
The diffractive contribution corresponds to dadiff/dMx — constant, while the non diffractive part (at 



large Mx) leads to a shape in 
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4 Measurement of the inclusive difFractive cross section HERA 



4.1 Inclusive diffraction as a leading twist process 

From observation of diffractive events, using the different techniques exposed above, the inclusive diffrac- 
tive cross section has been measured at HERA by HI and ZEUS experiments over a wide kinematic 
range [Il[2l[3lll]. Similarly to inclusive DIS, cross section measurements for the reaction ep — > eXp are 
conventionally expressed in terms of the reduced diffractive cross section, a^''^\ which is related to the 
measured cross section by 

a^^'\P,Q^xjp) . (2) 

At moderate inelasticities y, a^^'^' corresponds to the diffractive structure function ^2^''^'' to good ap- 
proximation. 

Fig. [7] illustrates a first result for the diffractive cross section as a function of W for different and 
Mx values. We notice that the diffractive cross section, ep —>■ epX, shows a hard dependence in the 
center-of-mass energy of the 7*p system W. Namely, we measure a W dependence of the form ~ W^'^ 
for the diffractive cross section, compatible with the dependence expected for a hard process. 

This first observation is fundamental and allows further studies of the diffractive process in the 
context of perturbative QCD (see next sections). The experimental selection of diffractive events is 
already a challenge but the discovery that these events build a hard scattering process is a surprise and 
makes the strong impact of HERA data into the field. Indeed, the extent to which diffraction, even in 
the presence of a hard scale, is a hard process, was rather unclear before HERA data. This has changed 
since then, with the arrival of accurate HERA data on diffraction in ep scattering and the realization 
that diffraction (measured to be a hard process) in DIS can be described in close analogy with inclusive 

DIS El Em. 

This is also confirmed in Fig. [HI where the ratio of diffractive to DIS cross sections is shown. This 
ratio is found to depend weakly on the Bjorken variable xbj (or W) at fixed values of the photon 
virtuality Q^. Thus, we can conclude that diffraction in DIS is a leading twist effect with logarithmic 
scaling violation in Q^, as for standard DIS. We discuss these results much further in the next sections. 

4.2 Recent results on inclusive diffraction at HERA 

Extensive measurements of diffractive DIS cross sections have been made by both the ZEUS and HI 
collaborations at HERA, using different experimental techniques [21 El II]- Of course, the comparison of 
these techniques provides a rich source of information to get a better understanding of their respective 
experimental gains and prejudices. 

In Fig. [9l the basis of the last ZEUS experimental analysis is summarized [3J. Data are compared to 
Monte-Carlo (MC) expectations for typical variables. The MC is based on specific models for signal and 
backgrounds, and the good agreement with data is proof that the main ingredients of the experimental 
analysis are under control. These last sets of data (Fig. ^ |3] contain five to seven times more statistics 
than in preceding publications of diffractive cross sections, and thus opens the way to new developments 
in data/models comparisons. 

A first relative control of the data samples is shown in Fig. [TOl where the ratio of the diffractive 
cross sections is displayed, as obtained with the LPS and the LRG experimental techniques. The mean 
value of the ratio of 0.86 indicates that the LRG sample contains about 24% of proton-dissociation 
background, which is not present in the LPS sample. This background corresponds to events like 
ep — >■ eXY, where y is a low-mass excited state of the proton (with My < 2.3 GeV). It is obviously 
not present in the LPS analysis which can select specifically a proton in the final state. This is the 
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Figure 7: The cross section of the diffractive process 7*p p'X, differential in the mass of the diffrac- 
tively produced hadronic system X (Mx), is presented as a function of the center-of-mass energy of 
the 7*p system W. Measurements at different values of the virtuality of the exchanged photon are 
displayed. 
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Figure 8: Left: Ratio of the diffractive versus the inclusive cross sections as a function of W for different 
values of and the diffractive mass Mx, derived from early ZEUS data. Right: Ratio of the diffractive 
versus total cross sections, as a function of x^j, derived from HI data for different values of and (3. 
A constant ratio of about 0.02 (2%) is observed for each bin of measurements. If we add up the five 
bins in (3 (for the bulk of the domain), we find immediately the average number of 10%. It gives the 
fraction of diffractive events on the total DIS sample (see text). 
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Figure 9: Left: Comparison of the distributions of data (dots) to those obtained with the Monte-Carlo 
(histograms) for typical variables in the LRG analysis. Right: Comparison of the distributions of data 
(dots) to those obtained with the Monte-Carlo (histograms) for typical variables in the LPS analysis. 

main background in the LRG analysis. Due to a lack of knowledge of this background, it causes a large 
normalization uncertainty of 10 to 15 % for the cross sections extracted from the LRG analysis. 

We can then compare the results obtained by the HI and ZEUS experiments for diffractive cross 
sections (in Fig. [H]), using the LRG method. A good compatibility of both data sets is observed, after 
rescaling the ZEUS points by a global factor of 13%. This factor is compatible with the normalization 
uncertainty described above. 

We can also compare the results obtained by the HI and ZEUS experiments (in Fig. [TTj) . using the 
tagged proton method (LPS for ZEUS and EPS for HI). In this case, there is no proton dissociation 
background and the diffractive sample is expected to be clean. It gives a good reference to compare both 
experiments. A global normalization difference of about 10% can be observed in Fig. [TTl which can 
be studied with more data. It remains compatible with the normalization uncertainty for this tagged 
proton sample. It is interesting to note that the ZEUS measurements are globally above the HI data 
by about 10% for both techniques, tagged proton or LRG. 

In Fig. [121 we compare the results using the LRG and the Mx methods, for ZEUS data alone. 
Both sets are in good agreement, which shows that there is no strong bias between these experimental 
techniques. The important message at this level is not only the observation of differences as illustrated 
in Fig. [m but the opportunity opened with the large statistics provided by the ZEUS measurements. 
Understanding discrepancies between data sets is part of the experimental challenge of the next months. 
It certainly needs analysis of new data sets from the HERA experiments. However, already at the present 
level, much can be done with existing data for the understanding of diffraction at HERA. 



4.3 Summary of recent results in one plot 

A summary of present measurements using LRG event selection is shown in Fig. [121 The ZEUS LRG 
data are extracted at the HI (3 and xjp values, but at different values. In order to match the 
Mn < 1.6 GeV range of the HI data, a global factor of 0.91 ±0.07, estimated with Pythia, is applied 
to the ZEUS LRG data in place of the correction to an elastic proton cross section. After this procedure, 
the ZEUS data remain higher than those of HI by 13% on average, as discussed above. The results of 
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Figure 10: Ratio of the diffractive cross sections, as obtained with the LPS and the LRG experimental 
techniques. The hues indicate the average value of the ratio, which is about 0.86. It implies that 
the LRG sample contains about 24% of proton dissociation events, corresponding to processes like 
ep eXY, where My < 2.3 GeV. This fraction is approximately the same for HI data (of course in 
the same My range). 
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Figure 11: Left: The diffractive cross sections obtained with the LRG method by the HI and ZEUS 
experiments. The ZEUS values have been rescaled (down) by a global factor of 13 %. This value 
is compatible with the normalization uncertainty of this sample. Right: The diffractive cross section 
obtained with the FPS (or LPS) method by the HI and ZEUS experiments, where the proton is tagged. 
The ZEUS measurements are above HI by a global factor of about 10%. 
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Figure 12: The diffractive cross sections obtained with the LRG method (full dots) compared with the 
results obtained with the Mx method (open symbols: FPC I and FPC II). All values are converted to 
My = Mp. 



the QCD fit to HI LRG data [2] is also shown (see next section). 



5 Diffraction and the resolved Pomeron model 
5.1 Theoretical considerations 

Several theoretical formulations have been proposed to describe the diffractive exchange. The purpose 
is to describe the blob displayed in Fig. [T]in a quantitative way, leading to a proper description of data 
shown in Fig. [71 

Among the most popular models, the one based on a point-like structure of the Pomeron assumes 
that the exchanged object, the Pomeron, is a color-singlet quasi-particle whose structure is probed in 
the reaction [5], E]- In this approach, diffractive parton distribution functions (diffractive PDFs) are 
derived from the diffractive DIS cross sections in the same way as standard PDFs are extracted from 
DIS measurements. It assumes also that a certain flux of Pomeron is emitted off the proton, depending 
on the variable xp, the fraction of the longitudinal momentum of the proton lost during the interaction. 
The partonic structure of the Pomeron is probed during the diffractive exchange O E] . 

In Fig. [3], we illustrate this factorization property and remind the notations for the kinematic 
variables used in this paper, as the virtuality of the exchanged photon, the center-of-mass energy of 
the 7*p system W and Mx the mass of the diffractively produced hadronic system X. It follows that 
the Bjorken variable xbj verifles xbj — Q'^/W'^ in the low xbj kinematic domain of the HI and ZEUS 
measurements {xsj < 0.01). Also, the Lorentz invariant variable /? deflned in Fig. [T]is equal to xbj/xp 
and can be interpreted as the fraction of longitudinal momentum of the struck parton in the (resolved) 
Pomeron. 

Because the short-distance cross section (7* — q) of hard diffractive DIS is identical to inclusive 
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Figure 13: Comparison between the HI and ZEUS LRG measurements after correcting both data sets to 
Mn < 1.6 GeV and applying a further scale factor of 0.87 (corresponding to the average normalization 
difference) to the ZEUS data. The measurements are compared with the results of the QCD fit prediction 
(see text). Further HI data at Xip = 0.03 are not shown. 
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DIS, the evolution of the diffractive parton distributions follows the same equations as ordinary parton 
distributions. Quantitatively, QCD factorization is expected to hold for F-f [S], E] and it may then be 
decomposed into diffractive parton distributions, ff, in a way similar to the inclusive F2, 

d^t d^t ^2'^^^'^'^)' 

where F2^i is the universal structure function for DIS on parton i, fx is the factorization scale at which 
//^ are probed and z is the fraction of momentum of the proton carried by the parton i. 

The QCD evolution equation applies in the same way as for the inclusive case. Fig. [8] is a simple 
experimental proof of this statement. For a fixed value of xp, the evolution in x and is equivalent 
to the evolution in P and Q^. 

If, following Ingelman and Schlein [6], one further assumes the validity of Regge factorization, 
may be decomposed into a universal Pomeron flux and the structure function of the Pomeron, 

^^^""^^plr''^ = /^/,(xp,t)Ff (AQ^) , (4) 

where the normalization of either of the two components is arbitrary. It implies that the xp and t 
dependence of the diffractive cross section is universal, independent of and P, and given by [2], O [6] 

/ , N 2ajp(0)-l 

In this approach, the mechanism for producing LRG is assumed to be present at some scale and 
the evolution formalism allows to probe the underlying partonic structure. The latter depends on the 
coupling of quarks and gluons to the Pomeron. It follows that the characteristics of diffraction are 
entirely contained in the input distributions at a given scale. It is therefore interesting to model these 
distributions. 



5.2 Diffractive parton densities 

In Fig. [THwe present the result for diffractive PDFs (quark singlet and gluon densities), obtained using 
the most recent inclusive diffractive cross sections presented in Ref. [2]. For each experiment (HI and 
ZEUS), we include measurements derived from Large Rapidity Gap (LRG) events in the QCD analysis. 

We follow the procedure described in Ref. [7], with previous ZEUS data. Note also that in all QCD 
fits, we let the global relative normalization of the data set as a free parameter (with respect to HI 
LRG sample) [7|. The typical uncertainties for the diffractive PDFs in Fig. ranges from 5% to 10% 
for the singlet density and from 10% to 25% for the gluon distribution, with 25% at large z (which 
corresponds to large /3 for quarks) [7j. Similar results have been obtained by the HI collaboration [2]. 

In order to analyze in more detail the large z behavior of the gluon distribution zG{z, = Ql) and 
give a quantitative estimate of the systematic error related to our parameterizations, we consider the 
possibility to change the gluon parameterization by a multiphcative factor (1 — zY (see Ref. [7]). If we 
include this multiplicative factor (1 — zY in the QCD analysis, we derive a value of u = 0.0 ± 0.5 (using 
the most recent data). Thus, we have to consider variations of u in the interval ±0.5 in order to allow 
for the still large uncertainty of the gluon distribution (mainly at large z values). The understanding 
of the large z behavior is of essential interest for any predictions at the Tevatron or LHC in central 
dijets production (see below). In particular, a proper determination of the uncertainty in this domain 
of momentum is necessary and the method we propose in Ref. [7] is a quantitative estimate, that can 
be propagated easily to other measurements. 

Of course, several checks need to be done to analyze the stability of the QCD fits procedure [7J. We 
present two of them below: 
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Figure 14: Singlet and gluon distributions of the Pomeron as a function of the fractional momentum 
of the Pomeron carried by the struck parton, derived from QCD fits on HI and ZEUS inclusive diffractive 
data (LRG)[2]. The parton densities are normalised to represent xp times the true parton densities 
multiplied by the flux factor at xp = 0.003 [7]. A good agreement is observed between both diffractive 
PDFs, which indicates that the underlying QCD dynamics derived in both experiments is similar. 




Figure 15: Left: Singlet and gluon distributions of the Pomeron as a function of the fractional 
momentum of the Pomeron carried by the struck parton, derived from QCD fits on HI data. Results 
are presented with = 3 GeV^ (full lines) and Q\ = 1.75 GeV^ (dashed lines), normalization follows 
the convention explained in Fig. [HI Right: Singlet and gluon distributions of the Pomeron as a 
function of z derived from QCD fits on HI data. Results are presented with Q^j„ = 4.5 GeV^ (full 
lines), Qfnin = 8.5 GeV^ (dashed lines) and Q^j^ = 12 GeV^ (dotted lines). 
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Figure 16: Predictions for x-pFj^ and = pD^pD as a function of /5 at = 30 GeV^ and xp = 10 ^ 

[7]. The dashed hne prediction refers to the diffractive PDFs analysis discussed in this part. Note that 
the longitudinal structure function Fp is directly related to the the reduced diffractive cross section: 
cr^*^^) = Fg^*-'^^ — .2 Fp^^\ Other curves represent dipole model calculations (see next sections). 



• We have checked the dependence of the diffractive PDFs on variations of the starting scale Qq in 
Fig. [15] (left). Very small changes are observed while changing the starting scale form 3 to 1.75 
GeV2. 

• We have checked the fit stability by changing the cut on Q^iny ^^e lowest value of of data to 
be included in the fit. The results are given in Fig. [15] (right), where we show the results of the 
fits after applying a cut on of 4.5, 8.5 and 12 GeV^. Differences are noticeable at small /? 
but well within the fit uncertainties. No systematic behavior is observed within Q^^^ variations. 

Then, an important conclusion is the prediction for the longitudinal diffractive structure function. 
In Fig. [16] [7j and [13 [8j, we display this function with respect to its dependence in /3 (Fig. [16] (a)) 
and the ratio R of the longitudinal to the transverse components of the diffractive structure function 
(Fig. [16] (b)). A comment is in order about the large (3 behavior. x^Fj^ is essentially zero at large (3 
from the pure QCD fits analysis. In fact, as illustrated in Fig. [16] and [TT] a non-zero contribution to 
the longitudinal structure function at large (3 corresponds to a twist-4, and is simply incorporated in 
a dipole model formulation of diffraction (see next sections). Here, we give the qualitative feature of 
this effect on the predictions for x^F^ . There is a significant difference between predictions with or 
without this twist-4 component in the region of large (3. However, the difference is negligible at low 
and medium /?, where the measurements are possible. Indeed, a first measurements has been realized 
by the HI collaboration [3], which is displayed in Fig. [18] It fits perfectly with the QCD fit prediction, 
as well as with the dipole prediction in the kinematic range accessible experimentally {(3 < 0.2). This 
Fig. [18] gives also directly the ratio of Fp versus F^. 
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Figure 17: Predictions for for xp = 10 ^ and = 10 GeV^ from fits with twist-4, dominant only 
at large (3, to the HERA data. The solid lines show predictions from pure (twist-2) QCD fits 
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Figure 18: Measurements of the diffractive longitudinal structure function multiplied by xp. The 
data are compared to QCD fits predictions (see text). Also shown for comparison is the value of F.^' as 
a dashed line. 
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6 Diffraction at the Tevatron and prospects for LHC 

6.1 Basics of diffraction at the Tevatron 

Once the gluon and quark densities in the Pomeron are known, it is easy to make predictions for the 
Tevatron (or the LHC) if one assumes that the same mechanism is the origin of diffraction in both 
cases. We assume the same structure of the Pomeron at HERA and the Tevatron and we compute 
as an example the jet production in single diffraction or double Pomeron exchange using the parton 
densities in the Pomeron measured at HERA. The interesting point is to see if this simple argument 
works or not, or if the factorization property between HERA and the Tevatron — using the same parton 
distribution functions — holds or not [9l [TOl [H] . In other words, we need to know if it is possible to use 
the parton distributions in the Pomeron obtained at HERA to make predictions at the Tevatron, and 
also further constrain the parton distribution functions in the Pomeron since the reach in the diffractive 
kinematical plane at the Tevatron and HERA is different. 

Theoretically, factorization is not expected to hold between the Tevatron and HERA ^ due to 
additional pp or pp interactions. For instance, some soft gluon exchanges between protons can occur 
at a longer time scale than the hard interaction and destroy the rapidity gap or the proton does not 
remain intact after interaction. The factorization break-up is confirmed by comparing the percentage of 
diffractive events at HERA and the Tevatron (10% at HERA and about 1% of single diffractive events 
at the Tevatron) showing already that factorization does not hold. This introduces the concept of gap 
survival probability, the probability that there is no soft additional interaction or that the event remains 
diffractive. 

The first experimental test of factorization concerns CDF data only. Fig. [19] shows the percentage 
of diffractive events as a function of x for different ^ bins and shows the same x-dependence within 
systematic and statistical uncertainties in all ^ bins supporting the fact that CDF data are consistent 
with factorization The CDF collaboration also studied the x dependence for different bins which 
leads to the same conclusions. 

A second step is to check whether factorization holds or not between Tevatron and HERA data. 
The measurement of the diffractive structure function is possible directly at the Tevatron. The CDF 
collaboration measured the ratio of dijet events in single diffractive and non diffractive events, which is 
directly proportional to the ratio of the diffractive to the standard proton structure functions, where 

^""^ Ratef^^ix)^ Fj^^ixY 

The comparison between the CDF measurement (black points, with systematics errors) and the expec- 
tation from the diffractive QCD fits on HERA data in full line is shown in Fig. [20] [TO] . 

We notice a discrepancy of a factor 8 to 10 between the data and the predictions from the QCD fit, 
showing that factorization does not hold. However, the difference is compatible within systematic and 
statistical uncertainties with a constant on a large part of the kinematical plane in /3, showing that the 
survival probability does not seem to be /3-dependent within experimental uncertainties. 

It would be interesting to make these studies again in a wider kinematical domain both at the 
Tevatron and at the LHC. The understanding of the survival probability and its dependence on the 
kinematic variables is important to make precise predictions on inclusive diffraction at the LHC. 

6.2 Discussion on the factorization breaking HERA/Tevatron 

In fact, from a fundamental point of view, it is natural that diffractive hard-scattering factorization does 
not apply to hadron-hadron collisions. Attempts to establish corresponding factorization theorems fail, 
because of interactions between spectator partons of the colliding hadrons. The contribution of these 
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interactions to the cross section does not decrease with the hard scale. Since they are not associated 
with the hard-scattering subprocess, we no longer have factorization into a parton-level cross section 
and the parton densities of one of the colliding hadrons. These interactions are generally soft, and we 
have at present to rely on phenomenological models to quantify their effects. 

The yield of diffractive events in hadron-hadron collisions is then lowered precisely because of these 
soft interactions between spectator partons (often referred to as re-interactions or multiple scatterings). 
They can produce additional final-state particles which fill the would-be rapidity gap (hence the often- 
used term rapidity gap survival). When such additional particles are produced, a very fast proton can 
no longer appear in the final state because of energy conservation. 

Diffractive factorization breaking is thus intimately related to multiple scattering in hadron-hadron 
collisions. We can also remark simply that the collision partners, in hadron-hadron reactions, are both 
composite systems of large transverse size, and it is not too surprising that multiple interactions between 
their constituents can be substantial. 

In contrast, the virtual photon in 7*p collisions has small transverse size, which disfavors multiple 
interactions and enables diffractive factorization to hold. According to our discussion, we may expect 
that for decreasing virtuality the photon behaves more and more like a hadron, and diffractive 
factorization may again be broken. 

6.3 Restoring factorization at the Tevatron 

The other interesting measurement which can be also performed at the Tevatron is the test of factoriza- 
tion between single diffraction and double Pomeron exchange. The results from the CDF collaboration 
are shown in Fig. [21] pDj. The left plot shows the definition of the two ratios while the right figure 
shows the comparison between the ratio of double Pomeron exchange to single diffraction and the QCD 
predictions using HERA data in full line. 

Whereas factorization was not true for the ratio of single diffraction to non diffractive events, fac- 
torization holds for the ratio of double Pomeron exchange to single diffraction. In other words, the 
price to pay for one gap is the same as the price to pay for two gaps. The survival probability, i.e. 
the probability not to emit an additional soft gluon after the hard interaction needs to be applied only 
once to require the existence of a diffractive event, but should not be applied again for double Pomeron 
exchange. 

To summarize, factorization does not hold between HERA and Tevatron as expected because of the 
long term additional soft exchanges with respect to the the hard interaction. However, experimentally, 
factorization holds with CDF data themselves and also between single diffraction and double Pomeron 
exchange which means that the soft exchanges do not depend on the hard scattering, which is somehow 
natural. 

6.4 Interest of exclusive events 

Once established some basics of the diffraction at the Tevatron, a fundamental topic concerns the 
analysis of exclusive events. A schematic view of non diffractive, inclusive double Pomeron exchange, 
exclusive diffractive events at the Tevatron or the LHC is displayed in Fig. [221 The upper left plot 
shows the standard non diffractive events where the Higgs boson, the dijet or diphotons are produced 
directly by a coupling to the proton and shows proton remnants. The bottom plot displays the standard 
diffractive double Pomeron exchange where the protons remain intact after interaction and the total 
available energy is used to produce the heavy object (Higgs boson, dijets, diphotons...) and the Pomeron 
remnants. We have so far only discussed this kind of events and their diffractive production using the 
parton densities measured at HERA. 
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Figure 19: Test of factorization within CDF data alone. The percentage of diffractive events is presented 
as a function of x for different ^ bins. The same x-dependence is observed within systematic and statis- 
tical uncertainties in all ^ bins, supporting the fact that CDF data are consistent with factorization [9]. 
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Figure 20: Comparison between the CDF measurements (Q^ = 75 GeV^, 0.035 < ^ < 0.095 and 
\t\ < 1 GeV^) of diffractive structure function (black points) with the expectation of the HERA (using 
first HI diffractive data) diffractive PDFs [7]. The large discrepancy both in shape and normalization 
between HERA predictions and CDF data illustrates the breaking of factorization at the Tevatron. 
Using the most recent measurements in QCD fits (and diffractive PDFs extraction) does not change 
this conclusion. 
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Figure 21: Restoration of factorization for the ratio of double Pomeron exchange to single diffractive 
events (CDF Coll.). Whereas factorization was not true for the ratio of single diffraction to non 
diffractive events, factorization holds for the ratio of double Pomeron exchange to single diffraction. 
The price to pay for one gap is the same as the price to pay for two gaps. 



There may be a third class of processes displayed in the upper right figure, namely the exclusive 
diffractive production. In this kind of events, the full energy is used to produce the heavy object 
(Higgs boson, dijets, diphotons...) and no energy is lost in Pomeron remnants. There is an important 
kinematical consequence. The mass of the produced object can be computed using roman pot detectors 
and tagged protons 

M = ^/M^. (7) 

We see immediately the advantage of these processes. We can benefit from the good roman pot resolution 
on and ^2 to get a good resolution on mass. It is then possible to measure the mass and the kinematical 
properties of the produced object and use this information to increase the signal over background ratio 
by reducing the mass window of measurement. It is thus important to know if this kind of events exist 
or not. 

In the following, we give some details of the search for exclusive events in the different channels 
which are performed by the CDF and DO collaborations at the Tevatron. Prospects for the LHC are 
then outlined. 

6.5 Search for exclusive events in Xc production 

For example, one way to look for exclusive events at the Tevatron is to search for the diffractive 
exclusive production of light particles like the x mesons. This would give rise to high enough cross 
sections - contrary to the diffractive exclusive production of heavy mass objects such as Higgs bosons 
— to check the dynamical mechanisms and the existence of exclusive events. Exclusive production 
of Xc has been studied by the CDF collaboration [TT] with an upper limit for the cross section of 
o'excipp —>■ p + J/ip + •y + p) ~49± 18{stat) ± 39{sys) pb, where the Xc decays into J/'if and 7, the 
J/\E' decaying itself into two muons. The experimental signature is thus two muons in the final state 
and an isolated photon, which is a very clear signal. 

Unfortunately, the cosmics contamination is difficult to compute and this is why the CDF collabora- 
tion only quotes an upper limit on the Xc production cross section. To know if the production is really 
exclusive, it is important to study the tail of inclusive diffraction which is a direct contamination of the 
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Figure 22: Scheme of non diffractive, inclusive double Pomeron exchange, exclusive diffractive events 
at the Tevatron or the LHC 



exclusive signal. The tail of inclusive diffraction corresponds to events which show very little energy in 
the forward direction, or in other words where the Pomeron remnants carry very little energy. This is 
why these events can be called quasi-exclusive. 

In Ref. [12], it is shown that the contamination of inclusive events into the signal region depends 
strongly on the assumptions on the gluon distribution in the Pomeron at high (3, which is poorly known 
as we mentioned in a previous section. Therefore, this channel is not conclusive concerning the existence 
of exclusive events. In the same spirit, the CDF collaboration also looked for the exclusive production 
of dilepton and diphoton [T3j. 

6.6 Search for exclusive events using the dijet mass fraction 

Another very important aspect of diffraction at the Tevatron is related to the diffractive production of 
dijet events in double Pomeron exchange (see Fig. [221) . The CDF collaboration measured the so-called 
dijet mass fraction in dijet events — the ratio of the mass carried by the two jets produced in the 
event divided by the total diffractive mass — when the antiproton is tagged in the roman pot detectors 
and when there is a rapidity gap on the proton side to ensure that the event corresponds to a double 
Pomeron exchange. 

The CDF collaboration has measured this quantity for different jet pt cuts [H]. In Fig. [231 we 
compare this measurement to the expectation coming from the structure of the Pomeron coming from 
HERA. For this sake, one takes the gluon and quark densities in the Pomeron measured at HERA as 
described in Ref. [7] and the factorization breaking between HERA and the Tevatron is assumed to 
come only through the gap survival probability (0.1 at the Tevatron). The comparison between the 
CDF data for a jet px cut of 10 GeV as an example and the predictions from inclusive diffraction is 
given in Fig. [231 l^ft. 

We also display in the same figure the effects of changing the gluon density at high j3 (by changing 
the value of the u parameter) and we note that inclusive diffraction is not able to describe the CDF 
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Figure 23: Left: Dijet mass fraction measured by the CDF collaboration compared to the prediction 
from inclusive diffraction based on the parton densities in the Pomeron measured at HERA. The gluon 
density in the Pomeron at high f3 was modified by varying the parameter u. Right: Dijet mass fraction 
measured by the CDF collaboration compared to the prediction adding the contributions from inclusive 
and exclusive diffraction. 



data at high dijet mass fraction, even after increasing the gluon density in the Pomeron at high (3 
(multiplying it by 1/(1 — j3)), where exclusive events are expected to appear [7| 115]. 

The conclusion remains unchanged when jets with > 25 GeV are considered [U [12]. Adding 
exclusive events to the distribution of the dijet mass fraction leads to a good description of data [H [T3] 
as shown in Fig. [231 right, where we superimpose the predictions from inclusive and exclusive diffraction. 

This study does not prove explicitly that exclusive events exist but shows that some additional 
component with respect to inclusive diffraction is needed to explain CDF data. Adding exclusive 
diffraction allows to explain the CDF measurement. To be sure of the existence of exclusive events, the 
observation will have to be done in different channels and the different cross sections to be compared 
with theoretical expectations [§. 

6.7 Prospects for LHC 

The search for exclusive events at the LHC can be performed in the same channels as the ones used at 
the Tevatron. Let us recall that a strong motivation for this idea is that heavy objects, like Higgs boson, 
could be produced in double pomeron exchange at the LHC [16j. In addition, some other possibilities 
benefiting from the high luminosity of the LHC appear. One of the cleanest ways to show the existence 
of exclusive events would be to measure the dilepton and diphoton cross section ratios as a function of 
the dilepton/diphoton mass [TB]. If exclusive events exist, this distribution should show a bump towards 
high values of the dilepton/diphoton mass since it is possible to produce exclusively diphotons but not 
dileptons at leading order as we mentioned in the previous paragraph. 

The motivation to install forward detectors at in ATLAS and CMS is then quite clear. In addition, 
it extends nicely the project of measuring the total cross sections in ATLAS and TOTEM by measuring 
hard diffraction at high luminosity at the LHC. Of course, this is a very challenging technical project. 

^ In Ref. [T3], the CDF data were also compared to the soft color interaction models. While the need for exclusive events 
is less obvious for this model, especially at high jet px, the jet rapidity distribution measured by the CDF collaboration is 
badly reproduced. This is due to the fact that, in the SCI model, there is a large difference between requesting an intact 
proton in the final state and a rapidity gap. 
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missing mass [GeV] 

Figure 24: Roman pot detector acceptance as a function of missing mass assuming a lOo" operating 
positions, a dead edge for the detector of 50 /im and a thin window of 200 /im. 



Without entering into details, a few technical issues can be discussed simply. Two locations for the 
forward detectors are considered at 220 and 420m respectively to ensure a good coverage in ^ or in 
mass of the diffractively produced object [16]. Installing forward detectors at 420m is quite challenging 
since the detectors will be located in the cold region of the LHC and the cryostat has to be modified to 
accommodate the detectors. In addition, the space available is quite small and some special mechanism 
called movable beam pipe are used to move the detectors close to the beam when the beam is stable 
enough. The situation at 220m is easier since it is located in the warm region of the LHC and both 
roman pot and movable beam pipe technics can be used. The AFP (ATLAS Forward Physics) project 
is under discussion in the ATLAS collaboration and includes both 220 and 420m detectors on both sides 
of the main ATLAS detector [TB] . 

To conclude on the diffraction at the LHC, the missing mass acceptance is given in Fig. [211 The 
missing mass acceptance using only the 220m pots starts at 135 GeV, but increases slowly as a function 
of missing mass. It is clear that one needs both detectors at 220 and 420m to obtain a good acceptance 
on a wide range of masses since most events are asymmetric (one tag at 220m and another one at 
420m). The precision on mass reconstruction using either two tags at 220m or one tag at 220m and 
another one at 420m is of the order of 2-4 % on the full mass range, whereas it goes down to 1% for 
symmetric 420m tags |16j . 
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Figure 25: Picture for the total cross section {'~f*p 7*p) in the dipole model. 

7 Diffraction and the dipole model 

7.1 Simple elements of theory 

The physical picture of hard diffraction at HERA is interesting in the proton rest frame and reminiscent 
of the aligned jet model. In the proton rest frame, at small xbj, the virtual photon splits into a qq 
pair long before it hits the proton [T71 UHl |19l [20] (see Fig. [25]). The qq wave-function of the virtual 
photon suppresses configurations in which one of the quarks carries almost all momentum. In fact, 
these configurations are the ones that give rise to a large diffractive cross section, just because the 
wave-function suppression is compensated by the large cross section for the scattering of a qq pair of 
hadronic transverse size off the proton. The harder of the two quarks is essentially a spectator to 
diffractive scattering. 

The scattering of the softer quark off the proton is non-perturbative and cannot be described by 
exchange of a finite number of gluons. Hence there is an unsuppressed probability that the softer 
quark leaves the proton intact. This explains simply the idea behind the leading twist nature of hard 
diffraction. The details of the scattering of the softer quark off the proton are encoded in the diffractive 
quark distribution. In a similar way, the qqg configuration in the virtual photon, in which the qq pair 
carries almost all momentum, gives rise to the diffractive gluon distribution. 

Also, in the simplest case, the colorless exchange responsible for the rapidity gap is modeled by the 
exchange of two gluons (projected onto the color singlet state) coupled to the proton with some form 
factor or to a heavy onium which serves as a model of the proton flEl \T9\ [20] . We focus the following 
discussion on dipole approaches of diffractive interactions, that follow exactly these ideas. 

Then, we can model the reaction in three different phases, as displayed in Fig. [26] -top-: 

(1) the transition of the virtual photon to the qq pair (the color dipole) at a large distance / ~ 
of about 10-100 fm for HERA kinematics, upstream the target, 

(2) the interaction of the color dipole with the target nucleon, and 

(3) the projection of the scattered qq onto the diffractive system X. 

7.2 Confrontation of HERA measurements to the dipole approach 

Following the arguments above, the inclusive diffractive cross section is described with three main 
contributions in dipole approaches. The first one describes the diffractive production of a qq pair from 
a transversely polarized photon, the second one the production of a diffractive qqg system, and the 
third one the production of a qq component from a longitudinally polarized photon (see Fig. [2B]-top-). 
In Fig. [23 -bottom-, we show that this approach, also called two-gluon exchange model gives a good 
description of the diffractive cross section measurements [TH] [191 [20] • 

One of the great advantage of the dipole model is that it provides a natural explanation of the 
rapidity gap formation. Another great advantage of the dipole formulation is that it provides a natural 
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Figure 26: Top: The qq and qqg components of the diffractive system. Bottom: The diffractive structure 
function xpF^^^^ is presented as a function of (3 for two values of Q^. The different components of the 
two-gluon exchange model are displayed (see text). They add up to give a good description of the data. 
The structure function xpF^^^^ is obtained directly from the measured diffractive cross section using 

the relation : ^'^^^"^^ °j^2 — ^^qT (1 — V + ^)F'^^^\xw>, x, Q"^), where y represents the inelasticity of the 
reaction. 



explanation of the experimental observation that cr'^^f-f/cr'^"^ ~ const as a function of energy W (see Fig. 
E]) [20] . Indeed, the dipole picture is valid in the frame in which the qq pair (dipole) carries most of the 
available rapidity Y ~ ln(l/x) of the system. The gluon radiation from the parent dipole can then be 
interpreted (in the large limit) as a collection of dipoles of different transverse sizes which interact 
with the proton. If the proton stays intact, diffractive events with large rapidity gap are formed. In such 
case, the diffractive system is given by the color dipoles and the diffractive exchange can be modeled 
by color singlet gluons exchange (two-gluon exchange) between the dipole and the proton (see Fig. [251 
and [26]-top-). When only the parent qq dipole forms a diffractive system, the diffractive cross section 
at t = reads 

= dhdz\^>\r,z,Q')\' a\x,r), (8) 

at I i=o Ibn J 

where is the well known light-cone wave function of the virtual photon, r is the dipole transverse 
size and 2 is a fraction of the photon momentum carried by the quark. Applying the qq dipole picture 
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Figure 27: The dipole cross section aqq in the saturation model as a function of dipole size r for different 
X (see text). 



to the total inclusive cross section, cr*°*, the following relation holds in the small- a; limit 



with the same dipole cross (3"(x, r) as in Eq. ([H]). This Eq. ^ is pictured in Fig. [23 
7.3 Saturation, concepts and practice 

In Eq. |8]and[9|, the parameterization of a{x,r) must be realized with caution [201 [19]. There are several 
features to consider. First, the density of gluons at given x increases with increasing Q^, as described in 
perturbative QCD. According to QCD evolution it also increases at given when x becomes smaller, 
so that the gluons become more and more densely packed. At some point, they start to overlap and thus 
re-interact and screen each other. Then, we enter a regime where the density of partons saturates and 
where the linear QCD evolution equations cease to be valid. To quantify these effects, a saturation scale 
Ql can be introduced, which also depends on x, such that for Q"^ ~ Ql{x) these effects of saturation 
become important. 

In practice, essential features of the saturation phenomenon are verified in the following parameter- 
ization for the dipole cross section first proposed in Ref. [20] 



where Qs{x) = Qo{x/xo)^^ is the saturation scale. In Fig. [271 we display the dipole cross section 
dependence of Eq. fllOp as a function of r at given x in this model. At small dipole size r ~ l/Q 
(large Q"^), the cross section rises following the relation a{x,r) oc r'^xg{x). At some value Rs{x) of r, 
the dipole cross section is so large that this relation ceases to be valid, and a{x, r) starts to deviate 
from the quadratic behavior in r. Therefore, Rs{x) = 1/Qs{x) represents a typical saturation scale. As 
r continues to increase, a{x, r) eventually saturates at a value typical of a meson-proton cross section. 
For smaller values of x, the initial growth of agg with r is stronger because the gluon distribution is 
larger. The target is thus more opaque and saturation sets in at lower r. 

Parameters of the dipole cross section of Eq. ([TOj) are obtained from the analysis of inclusive data, 
and then can be used to predict diffractive cross section in DIS, and even more processes as we discuss in 
the next section. An important aspect of Eq. (ITOl) . in which r and x are combined into one dimensionless 




(9) 



a{x,r) = (To{l - exp{-r'^Ql(x))} 



(10) 
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variable rQs{x), is what is called geometric scaling, a new scaling property in inclusive DIS at small x. 
In Ref. [20], it has been shown to be valid for the total cross section (see Fig. 17.31) . 
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Figure 28: The total cross section as a function of r = Q'^ /Ql for x < 0.01. The saturation 

regime is reached when Qs^Q- For inclusive events in deep inelastic scattering, this feature manifests 
itself (as displayed) via the geometric scaling property: instead of being a function of Q^/Ql and x 
separately, the total cross-section is only a function of t = Q'^/Q1{x), up to large values of r. 



It happens that diffraction in DIS is an ideal process to study parton saturation since this process is 
especially sensitive to the large dipole contribution, r > 1/Qs{x) [20j. Unlike inclusive DIS, the region 
below, r < 1/Qs{x), is suppressed by an additional power of This makes obviously diffractive 

interactions very important for tracting saturation effects. As already mentioned, the dipole cross 
section with saturation (see Eq. ffTOl) ) leads in a natural way to the constant ratio (up to logarithms) 



We can present very simply the main elements of the calculation that bring this result. Indeed, the 
photon wave function, in Eq. ([8]), favors small dipoles (small r ~ 1/Q), which gives 



\d\dz\^\r^zM')?o\x^r)^^\ ^a\x,r) 
dt \t=o IQn I ./i/o2 



On the other hand, the dipole cross section favors relatively large dipoles, with a{x, r) ~ r^. However, 
as discussed above in the building of Eq. (ITOl) . at sufficiently high energy, saturation cuts off the large 
dipoles already on the semi-hard scale 1/Qs- This leads to 

da'^^ff 1 f'/^Ur^ Ql{x) 



at \t=o Ji/Q2 
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Figure 29: The unified picture of Compton scattering, diffraction excitation of tlie plioton into liadronic 
continuum states and into tlie diffractive vector meson 

and it follows immediately that is a constant of x at fixed values of Q^. This result is illustrated 
experimentally at the beginning of this review, in Fig. [81 

With Eq. [TOl we have also introduced above an interesting consequence of the dipole model for the 
total cross section, the geometric scaling property. Namely, the total cross section does not depend on 
X and independently but can be expressed as a function of a single variable r = Q^/Ql{x) [20] . 
This property has also been shown recently to be verified under minimal assumptions for all diffractive 
processes [21] (see Fig. [30ll . The experimental confirmation of this relation is an interesting piece of 
evidence that saturation effects are (already) visible in the inclusive diffractive DIS data. Extensions 
of these ideas at non-zero t values, rooted on fundamental grounds, have also been recently derived 
|19j . This provides essential perspectives to understand the transverse degrees of freedom which are 
discussed in the next sections. 



7.4 Towards a common description of all diffractive processes 

Let us mention that one of the great interest of the dipole model in its two-gluon exchange formulation 
is that it provides a unified description of different processes measured in 'y*p collisions at HERA: 
inclusive 7*p X, diffractive 7*p ^ X p' and (diffractive) exclusive vector mesons (VM) production 
■j*P VM p' (see Fig. [29]). In the last case, the step (3) described in Fig. [25] consists in the 
recombination of the scattered pair qq onto a real VM (as J/\E', p°, 0,.--) [^ [22] [21] [25] or onto a 
real photon for the reaction 7*p — > 7 p'. This last process is called deeply virtual Compton scattering 
(DVCS) [26] [27]. Also, we understand immediately the fundamental interest of exclusive VM production 
to clarify the generic mechanism of diffractive DIS. Indeed, the scales involved in the VM process can 
act as triggers to isolate when the virtual photon fluctuates mainly into small-size (small r) qq pair 
configurations, or mainly into large-size configurations. For example, a small-size qq dipole is most 
likely to be produced if the virtual photon is polarized longitudinally or if the dipole is built with heavy 
quarks. Therefore, we can already state that exclusive J/\E' production is a good candidate for a hard 
diffractive process fully calculable in perturbative QCD. We discuss completely these ideas in the next 
section. 
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Figure 30: The diffractive cross section /5 daj-jj* ^/dp from HI and ZEUS measurements, as a function 
of Td = Q'^/Q'iixp) in bins of /3 for Q"^ values in the range [5; 90] GeV^ and for xp < 0.01 [21j. The 
geometric scahng is confirmed with a good precision. 
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8 Exclusive particle production at HERA 



8.1 Triggering the generic mechanism of difFractive production 

There is a long experimental and theoretical history to the study of vector meson production, revived 
with the advent of HERA. On the experimental side, the important result is that the cross sections for 
exclusive vector meson production rise strongly with energy (if a hard scale is present) when compared 
to fixed target experiments. A compilation of experimental measurements are shown in Fig. ED [221 
[23t [211 [25]. We observe some statements mentioned briefly at the very end the previous section. For 
example, for J/ip exclusive production, the W dependence of the cross section is typical of a hard 
process. Indeed, the mass of the J/ip plays the role of the large scale, which mainly triggers small- 
size (small r) qq pair configurations of the initial virtual photon, which then build the hard process. 
If we follow the discussion of the previous section, we can also easily write the above argument at a 
quantitative level. Indeed, VM cross sections in the (hard) perturbative regime, 7*p — > VM p', depend 
on the square of the gluon density in the proton. A first approximation of the cross section can then 
be written as 

da 
H 

where the dependence on the meson structure is in the parameter 



(7*iV VN) = An'Tymya^iQf-i^^i^^^^ , (13) 



Vv = ^ I -^^0^(^) ( / rf#^(^) ) (14) 



2 J z{l-z) 

and (p^iz) is the leading- twist light-cone wave function. 

Fig. [31] presents also interesting features that we can comment at this level of the discussion. It 
shows the transition from soft to hard processes, using the mass of the VM as a trigger. From the 
lightest one, p°, up to the T, Fig. [3T] shows ct(7p Vp) as a function of W. For comparison, the total 
photoproduction cross section, atoti'jp), is also shown. The data at high W can be parameterized as 
W^, and the value of 5 is displayed in Fig. [31] for each reaction. One sees clearly the transition from a 
shallow W dependence for low mass VM (soft) to a steeper one as the mass of the VM increases (hard) 

[221 [221 [21 Eg. 

An interesting phenomenon is observed for the DVCS cross section (see Fig. [32]) . which presents the 
same hard W dependence as for the J/-?/' [26] [2Z], with a (zero mass) photon in the final state. It does 
not seem to follow the logic of the above argument and we come back later of this point. Obviously, 
Fig. [nH displays only one aspect of the problem, using the mass of the VM as the scale trigger for the 
soft-hard diffractive process. It is clear that the scale is also particularly well suited, always for the 
exclusive electroproduction of light vector mesons [221 1221 1211 [22] and DVCS [23 [2Z] ■ The soft-hard 
transition can be observed experimentally in different ways when varying Q^: 

(1) In the change of the logarithmic derivative 6 of the process cross section a with respect to the 7*p 
center-of-mass energy W {a ^ W^). We expect a variation from a value of about 0.2 in the soft 
regime (low values) to 0.8 in the hard one (large values). 

(2) In the decrease of the exponential slope b of the differential cross section with respect to the 
squared- four- momentum transfer t {da /dt ~ e"'*'*'), from a value of about 10 GeV~^ to an asymp- 
totic value of about 5 GeV~^ when the virtuality of the photon increases. 

We illustrate this procedure on recent data on p° production [22]- The cross section a{'y*p — ^ p^p) 
is presented in Fig. [33] as a function of W, for different values of Q^. The cross section rises with W in 
all bins. The same conclusion holds for DVCS, as already discussed and shown in Fig. [321 ^26^ i27j . 
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Figure 31: dependence of the exclusive vector meson cross section in photoproduction, ai^jp — > Vp). 
The total photoproduction cross section is also shown. The lines are the fit result of the form to 
the high energy part of the data. 
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Figure 32: The DVCS cross section, a'^ p^ip ^ function of W for different values. The solid 
lines are the results of a fit of the form a'^'p^'^P oc . The values of 5 and their statistical uncertainties 
are given in the figure. 
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Figure 33: W dependence of the cross section for exclusive p° electroproduction, for different values, 
as indicated in the figure. The lines are the fit results of the form to data. 
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Figure 34: A compilation of values of 6 from fits of the form for exclusive VM electroproduction, 
as a function of + M^. It includes also the DVCS results. 
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A compilation of values of 6 from DVCS and VM measurements are presented in Fig [311 Results 
are plotted as a function of + M^, where M is the mass of the vector meson (equal to zero in case 
of DVCS). We observe a universal behavior, showing an increase of S as the scale becomes larger. The 
value of 6 at low scale is the one expected from the soft Pomeron intercept, while the one at large scale 
is in accordance with twice the logarithmic derivative of the gluon density with respect to W. 

8.2 Deeply virtual Compton scattering 




Figure 35: Diagrams illustrating the DVCS (left) and the Bethe-Heitler (middle and right) processes. 
The reaction studied receives contributions from both the DVCS process, whose origin lies in the strong 
interaction, and the purely electromagnetic Bethe-Heitler (BH) process, where the photon is emitted 
from the positron. The BH cross section can be precisely calculated in QED using elastic proton form 
factors. 

Let us comment in more details the analysis of the DVCS signal, which we discuss in a different 
context in further sections. The DVCS process, ep ep7, also receives a contribution from the purely 
electromagnetic Bethe-Heitler (BH) process, where the photon is emitted from the electron, as displayed 
in Fig. [351 

Let us notice that the final state for DVCS (QCD process) and BH (QED process) are identical. 
This means that both processes interfere, which is of fundamental interest in the next sections. In 
this part, we only use the fact that the BH cross section is precisely calculable in QED and can be 
subtracted from the total process rate to extract the DVCS cross section. Of course, only if the BH 
contribution is not dominating the process rate and if the (integrated) interference term is negligible. 
Otherwise, the subtraction procedure would be hopeless. It is the case at low xbj, and then for HI and 
ZEUS experiments, the DVCS contribution can be measured directly. Fig. [361 presents the different 
contribution (for the scattered electron variables), after the experimental analysis of the reaction ep —>■ 
ep7. 

We observe that DVCS and BH contributions are of similar size and thus, the BH contribution can 
be subtracted with a systematic uncertainty determined from a specific experimental study. In Fig. [371 
we present the DVCS cross sections, 7*p 7p, obtained over the full kinematic range of the analysis 
[261 [27], as a function of and W. The behavior in W has been discussed qualitatively above, it 
corresponds to the dependence characteristic for a hard process. The dependence, measured to be 
in ~ l/Q^ in Fig. ^71 is also understandable quahtatively. 

Indeed, following the discussion of the previous section (see Eq. ([T2l) ). we expect a behavior of the 
imaginary DVCS amplitude {'~f*p — > 7p) in 
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Figure 36: Distributions of the energy and polar angle of the scattered electron. The data are compared 
with Monte Carlo expectations for elastic DVCS, elastic and inelastic BH and inelastic DVCS (labelled 
DISS. p). All Monte Carlo simulations are normalized according to the luminosity of the data. The 
open histogram shows the total prediction and the shaded band its estimated uncertainty. 



which leads to a DVCS cross section of the form 

With this expression, we find again the qualitative behavior in W . Interestingly also, the measured 
dependence in ~ is smaller than expected from this relation. In fact, to describe qualitatively the 
observed DVCS cross section, we must consider a parameterization in 

after introducing a term in [Q'^]"' in the expression of the DVCS cross section. The term in [Q'^Y' 
reminiscent from the QCD evolution of the DVCS amplitude (QCD evolution of the gluon/sea distri- 
butions). The experimental observation in a ~ is compatible with 7 ~ 1/2 (using our notations). 
Of course, we do not stay at this qualitative understanding and we describe quantitative estimates of 
the DVCS cross sections in the following. 

A comment is in order concerning the W dependence of DVCS. It reaches the same value of 5 as 
in the hard process of J/ip electroproduction. Given the fact that the final state photon is real, and 
thus transversely polarized, the DVCS process is produced by transversely polarized virtual photons, 
assuming s-channel helicity conservation. The steep energy dependence thus indicates that the large 
configurations of the virtual transverse photon are suppressed and the reaction is dominated by small 
qq configurations (small dipoles), leading to the observed perturbative hard behavior. A similar effect 
is observed for p° production [22J . 



8.3 Saturation in exclusive processes 

Coming back to the discussion of the previous section about saturation, we can mention also that among 
diffractive interactions, exclusive vector meson production and DVCS are probably the best processes 
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Figure 37: The DVCS cross section as a function of at W = 82 GeV and as a function of W at 
= 8 GeV^. The inner error bars represent the statistical errors, the outer error bars the statistical 
and systematic errors added in quadrature. 

to study saturation effects in DIS since the transverse size of the qq pair forming a meson is controlled 
by the vector meson mass with < r >~ 1 / a/ My + Q'^ . Thus we expect saturation effects to be more 
important for larger (lighter) vector mesons. An interesting consequence of this feature is illustrated in 
Fig. [38|, where we show that VM and DVCS process exhibit the property of geometric scaling [21j. This 
illustrates that this qualitative discussion (related to Eq. f[T^ ) gives the main elements of understanding 
of the DVCS and VMs cross sections dependences. More generally, as for all other diffractive processes 
presented in this review, it means that the mechanism included in the parameterization of the dipole 
cross section of the form written in Eq. (fTOj) is correct and predictive. 

In recent works, it has been shown that dipole models can be extended at non-zero t values, with 
a refined definition of the saturation scale [19]. Then, the geometric scaling property is predicted to 
manifest itself in exclusive vector meson production and deeply virtual Compton scattering (DVCS), 
also at moderate non-zero momentum transfer. In Fig. [39l we compare data with predictions of Ref. 
[To] , We observe the very good agreement between data and predictions. 
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Figure 38: The p, J/\l/ and (f) production cross-sections (Tvm^^^ and the DVCS cross-section (t]jvcs^ 
from HI and ZEUS measurements, as a function of ry = {Q^ + My^j)/Ql{x]p) and for xp <0.01 [21]. 
Each process verifies the geometric scahng property. 




Figure 39: Fit resuhs for the p, and J/\I/ differential cross-section. 
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9 Nucleon tomography 



9.1 t dependence of exclusive difFractive processes revisited 

With t = {p — p'Yi the measurement of the VM and DVCS cross section, differential in t is one of 
the key measurement in exclusive processes. A parameterization in da/dt ~ e~*'*', as shown in Fig. 
HOI gives a very good description of measurements. In addition, in Fig. HOI we show that fits of the 
form da/dt ~ e"'''*' can describe DVCS measurements to a very good accuracy for different and W 
values. The same conclusions hold in the case of VM production. That's the reason why we use this 
parameterization of the t dependence, with a factorized exponential slope b, to describe the HERA data 
on DVCS or VM production at low xbj- Note that this parameterization 

Concerning the interpretation, we have already briefly mentioned the importance of the observation 
of the decrease of the exponential slope b, from a value of about 10 GeV~^ to an asymptotic value of 
about 5 GeV~^, when the virtuality of the photon increases (see Fig. HTI) . The resulting values of b 
as a function of the scale + are plotted in Fig. |4T1 




Figure 40: The DVCS cross section, differential in t, for three values of expressed atand for three 
values of W. The solid lines represent the results of fits of the form e"^'*'. 

A qualitative understanding of this behavior is simple. Indeed, b is essentially the sum of a component 
coming from the probe in 1/ a/QM-M^^ and a component related to the target nucleon. Then, at 
large or large My^j, the b values decrease to the solely target component. That's why in Fig. SH 
we observe that for large or for heavy VMs, like J/ip., b is reaching a universal value of about 5 
GeV~^, scaling with Q"^ asymptotically. This value is related to the size of the target probed during the 
interaction and we do not expect further decrease of b when increasing the scale, once a certain scale is 
reached. 

To understand this shape of more quantitatively, we need to define a function that generalizes 
the gluon density which appears in Eq. (fT3!) at non-zero t values. That's why, we define a generalised 
gluon distribution Fg which depends both on x and t (at given Q"^). From this function, we can 
compute a gluon density which also depends on a spatial degree of freedom, a transverse size (or impact 
parameter), labeled i?^, in the proton. Both functions are related by a Fourier transform 
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At this level of the discussion, there is no need to enter into further details concerning these functions. 
We just need to know that the functions introduced above define proper (generalized) PDFs, with gauge 
invariance and all the good theoretical properties of PDFs in terms of operator product expansion. In 
fact, they are rooted on fundamental grounds [28], that we develop in further sections (without heavy 
formalism) . 



9.2 Extracting the transverse distribution of the quarks and gluons 

From the Fourier transform relation above, the average impact parameter (squared), of the dis- 

tribution of gluons g{x, Rj_) is given by 



(4) - 



/ d^R^g{x,Rj 



J dm^g{x,R, 



d_ 
di 



LF,(x,0) 



26, 



(16) 



i=0 



where b is the exponential t-slope. In this expression, (r|.) is the transverse distance between the 
struck parton and the center of momentum of the proton. The latter is the average transverse position 
of the partons in the proton with weights given by the parton momentum fractions. At low xbj, the 
transverse distance defined as \/Jr^ corresponds also to the relative transverse distance between the 
interacting parton (gluon in the equation above) and the system defined by spectator partons. Therefore 
provides a natural estimate of the transverse extension of the gluons probed during the hard process. 

In other words, a Fourier transform of momentum to impact parameter space readily shows that the 
t-slope b is related to the typical transverse distance in the proton. This t-slope, b, corresponds exactly 
to the slope measured once the component of the probe itself contributing to b can be neglected, which 
means at high scale: or Myj^. Indeed, at high scale, the qq dipole is almost point-like, and the t 
dependence of the cross section is given by the transverse extension of the gluons in the proton for a 
given XBj range. 



9.3 Comments on the physical content of (r^) 

A short comment is in order concerning the fundamental relation (fT6!) for DVCS at HERA (at low xbj)- 
Does it make sense to keep only the gluon distribution in this expression or do we need to consider 
also sea quarks? This issue can be addressed simply by coming back to Eq. (IT^ . where we have 
approximated the imaginary DVCS amplitude (7*^ 7p) in 

)2 /_ ,„„ 'Za 



Let us give first a more general form to this formula, keeping the tracks of the photon wave functions 

ImA = j (frdz'^*{r,z,Ql = Q'^)'^{r,z,Ql = 0)a{x,r), (17) 

where '^*{r, z,Ql = Q"^) is the wave function for the virtual photon and \I'(r, 2;,Q2 = 0) the real 

photon. Also, following the previous discussion on the dipole cross section, we can write: a{x,r) ~ 
cror'^Qs{x,r)'^, with 

Qs{x,r) — Qoi—) , 

where Rp is the proton radius. We conclude immediately that the imaginary part of the DVCS amplitude 
is dependent on the gluon density convoluted by the photon (virtual and real) wave functions. It gives 
the rationale behind formula (fT6l) . 
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Of course, this is a matter of representation. In the Eq. f|T7|) . we write the photon-gluon interaction 
through a quark loop, with a virtual photon fluctuating in a qq pair, which is exactly the dipole qq 
component entering in \l/*(r, = Q^) (see also Fig. [22]). In other words, at low xbj {xej — 10~^), 
the idea is that quarks (sea quarks) are produced by gluons. 

Then, the dipole formalism, summarized in Eq. (fT7|) or Fig. [29], provides a very powerful expression 
of this behavior. Of course, in other formalisms, that we present latter, we can express the cross sections 
at the level of the photon-quark interaction and thus consider directly the sea quark distribution. 

9.4 Experimental results 

DVCS results lead to = 0.65 ± 0.02 fm at large scale > 8 GeV^ for xsj ^ 10"^ [26]. This value 
is smaller than the size of a single proton, and, in contrast to hadron-hadron scattering, it does not 
expand as energy W increases (see Fig. H2]) . Then, we can parametrize the measured b values displayed 
in Fig. [32] in the form of a Pomeron trajectory: b = + 2a' In ^— . We obtain that the a' value, which 
is characteristic of the energy dependence of the trajectory, is close to zero. 

This is not useless to recall that this observation is extremely challenging on the experimental 
analysis side. We are dealing with nano-barn cross sections, that we measure as a function of t, and 
finally, we measure the energy dependence of this behavior in t. Of course, the gain is important. In 
particular, the great interest of the DVCS is that the t dependence measured is free of effects that could 
come from VM wave functions (in case of VMs) and then spoil (to a certain limit) the interpretation of b 
described above. Thus, with DVCS, we have the advantage to work in a controlled environment (photon 
wave functions) where the generic Eq. (1161) can be applied to the measurement (almost directly) and 
must not be corrected with effects arising from VMs wave function. 

It is obviously very interesting to extend the result presented in Fig. [12] to all VMs. Indeed, we 
can study the W dependence of dcr/dt and extract the energy dependence as done above for all VMs, 
using b = bo + 2a' In Results are presented in Fig. [33 (bottom). Values are plotted as a function 

of + M^. We observe that the values of a' tends to decrease with the scale. In particular, the 
measurement of a' done for the J/\l/ [23], leading to a small value for a', is well compatible with the 
DVCS result [26]. 

A short comment can be done qualitatively on such small a' value. We can rephrase this observation 
as an evidence of no shrinkage of da/dt in the process 7*p J/'^p or 7*p — > 7p. Looking at the 
diagram describing two gluon exchange in Fig. [331 the virtual photon fluctuates into two high 
quarks. Although in the diagram there are only two gluons linked to the proton, we actually have a 
whole ladder due to the large rapidity range available at these high W energies (see Fig. [33]) • From the 
virtual photon vertex down to the proton, the average kx of the gluons gets smaller, the configuration 
larger and we enter the region of low kx physics governed by non-perturbative QCD. This process is 
called Gribov diffusion. Thus a process can start as a hard process at the photon vertex but once 
it couples to the proton it gets a soft component which makes the process non calculable in pQCD. 
The average kr of the partons in the process can be estimated by the slope of the trajectory since 
a' ~ 1/ < kx >■ 

The fact that no shrinkage is observed indicates that Gribov diffusion is not important in this process 
at the presently available W values, and the average kx remains large. Such a behavior is expected for 
hard processes, where a' <C 0.25 GeV~^. The experimental results for exclusive J/\& production and 
DVCS confirm that both processes are fully calculable in perturbative QCD. 

9.5 Link with LHC issues 

Let us finish this section by a comment making the link with LHC issues. Indeed, the correlation 
between the transverse distribution of partons and their momentum fraction is not only interesting 
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Figure 41: A compilation of the values of the slope b as a function of + for various exclusive 
processes including the present DVCS measurement. The inner error bars represent the statistical 
uncertainty while the outer error bars the statistical and systematic uncertainties added in quadrature. 
Note that the latest t slope measurement of DVCS by the ZEUS collaboration [27] is shown. It falls 
below (1 sigma effect) the HI measurement at a comparable value |26]. The main result does not 
change: at large Q^, exponential t slopes converge to a scaling value (see text) ans this is a common 
trend for all VM processes. However, at low [Q"^ ~ 3 GeV^), ZEUS result indicates the absence of 
effects in 6 ~ l/Q^ (from the probe) for DVCS, whereas HI result shows a behavior comparable to p 
at this Q^, with a clear influence of the probe to the building of the measured b. 
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Figure 42: The logarithmic slope of the t dependence for DVCS as a function of W. 
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Figure 43: Values of the intercept and slope of the effective Pomeron trajectory as a function of 
+ M^, as obtained from measurements of exclusive VM electroproduction. 
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Figure 44: A diagram describing a gluon ladder in a diffractive process. 
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Figure 45: a: Graph with a single hard interaction in a hadron-hadron collision. The impact parameters 
h\ and 62 are integrated over independently, b: Graph with a primary and a secondary interaction. 
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from the perspective of hadron structure, but also has practical consequences for high-energy hadron- 
hadron collisions. Consider the production of a high- mass system (a dijet or a heavy particle). For the 
inclusive production cross section, the distribution of the colliding partons in impact parameter is not 
important: only the parton distributions integrated over impact parameters are relevant according to 
standard hard-scattering factorization (see Fig. H5](a)). There can however be additional interactions in 
the same collision, especially at the high energies for the Tevatron or the LHC, as shown in Fig. HST b). 
Their effects cancel in sufficiently inclusive observables, but it does affect the event characteristics and 
can hence be quite relevant in practice. In this case, the impact parameter distribution of partons must 
be considered. 

The production of a heavy system requires large momentum fractions for the colliding partons. A 
narrow impact parameter distribution for these partons forces the collision to be more central, which 
in turn increases the probability for multiple parton collisions in the event (multiple interactions). 

10 Generalised parton distributions 
10.1 A brief introduction in simple terms 

We have already defined in a previous section a first form for a generalized gluon distribution. In 
this part, we move into further details and explain the wide experimental field opened in the area of 
generalized parton distributions. 

First, a short contrarian comment: DIS can not be considered as the continuation of the original 
Rutherford experiment. Indeed, Rutherford measured that the nucleus is concentrated in a very small 
part of the atom, and, as far as we consider only PDFs, we have no possibility to explore the spatial 
structure of the nucleon. The reason is that in the infinite momentum frame picture, the light-cone 
description of the Feynman parton model does not explore the space-time location of partons. In other 
words, within the infinite momentum frame description, the variable xbj has no direct relation to the 
space coordinate of a parton but is related to a combination of the energy and momentum of this parton. 

In the previous section, we have shown that data on exclusive particle production can give access to 
the spatial distribution of quarks and gluons in the proton at femto-meter scale. Then, we have defined 
functions, which model this property (for gluons) through the relation 



Of course, a similar relation holds for quarks, linking the two functions q{x, R±; Q"^) and Fg{x,t = 
— A^; Q"^). The general framework for this physics is encoded in the so-called generalized parton distri- 
butions (GPDs). 

We already know that the reconstruction of spatial images from scattering experiments by way of 
Fourier transform of the observed scattering pattern is a technique widely used in physics, for example, 
in X-rays scattering from crystals. In simple words, what we have done experimentally is that we 
have extended this technique to the spatial distribution of quarks and gluons within the proton, using 
processes that probe the proton at a tiny resolution scale. Of course, as already mentioned, working at 
a femto-meter scale with nano-barn cross sections is very challenging from the experimental front. We 
have achieved this and it immediately opens a way in the ambitious program of mapping out the GPDs. 
We come back below in a more systematic way on different aspects of that program that requires a large 
amount of experimental informations, for which future programs at JLab and CERN are appealing. 

Before coming back to the experimental side, we can present a short overview of GPDs, in simple 
terms. It is interesting, even for an experimentalist, as it clarifies the Fourier transform relation discussed 
above and makes more transparent the goals for the future. For complete reviews, see Ref. [28l [29| [30]. 
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Figure 46: Picture of a GPD and its variables. The momentum fractions x and ^ refer to the average 
hadron momentum \{p + p')- Note that x is an internal variable and is not equal to xbj- However, 
there is a relation between the skewing variable ^ and xbj, ^ = XBj/i^2 — xbj)- 



GPDs are defined through matrix elements {p'\0\p) between hadron states \p') and \p), with non-local 
operators O constructed from quark and gluon fields. From this expression, we understand why GPDs 
are directly related to the amplitude for VM or real gamma exclusive production. For unpolarized 
quarks there are two distributions H'^{x,C,,t) and E'^{x,^,t), where x and ^ are defined in Fig. HHl 
The former is diagonal in the proton helicity, whereas the latter describes proton helicity flip. For 
p = p' and equal proton helicities, we recover the diagonal matrix element parameterized by usual 
quark and antiquark densities, so that H'^{x,0,0) = q{x) and H'^{—x,0,0) = ~q{x) for x > 0. Note 
that the functions of type E are not accessible in standard DIS, as it corresponds to matrix elements 
{p' , s'\0\p, s) with s 7^ s'. Even in DVCS-like analysis, it is very difficult to get a sensitivity to these 
functions, as in most observables, their contributions are damped by kinematic factors of orders \t\/Mp, 
with an average \t\ value in general much smaller that 1 GeV^. Then, till stated otherwise, our next 
experimental discussions are concentrated on the determination of GPDs of type Hq or Hg. We come 
back later on this point and show specific cases where E-type functions can be accessed and why this 
is an important perspective. 

10.2 Fundamental relations between GPDs and form factors 

An interesting property of GPDs, which lightens their physics content, is that their lowest moments 
give the well-known Dirac and Pauli form factors 

[ dxH\x,^,t) = I dxE\x,i,t) = F2(t), (18) 

where denotes the fractional quark charge. It means that GPDs measure the contribution of 
quarks/gluons, with longitudinal momentum fraction x, to the corresponding form factor. In other 
words, GPDs are like mini-form factors that filter out quark with a longitudinal momentum fraction x 
in the proton. Therefore, in the same way as Fourier transform of a form factor gives the charge distri- 
bution in position space, Fourier transform of GPDs (with respect to variable t) contains information 
about the spatial distribution of partons in the proton. 

10.3 New insights into proton imaging 

This discussion clarifies also the Fourier transforms, that can relate g{x, R±; Q'^) and Fg{x, t = — A^^; Q^) 
or q{x, R±; Q'^) and Fq{x,t = —A']_;Q'^). We have already discussed these functions and from their 
relations, it follows that q{x, Rj_; Q^) is the probability density to find a quark with momentum fraction 
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a; at a transverse distance R± from the (transverse) center of momentum of the proton. More formal 
discussions can be found in Ref. [2B]. 

Exactly, what must be confronted with the proton radius is not \/r^ but — xbj), which 

does not change our result with xbj — 10^^ (v^^ ~ 0-65 ± 0.02 fm), but must be taken into account 
for fixed target kinematics at larger xbj- In particular, at very large xbj {xbj 1), the struck quark 
is carrying almost the entire proton momentum, thus its relative distance to the center of momentum 
of the proton obviously tends to zero. This means that \/r^ tends to zero (by definition). In order to 
keep finite the ratio -^/r^/(l — xbj), we can conclude that the asymptotic form of \/r^ at large xbj is 
likely in (1 — XBjY- 

The distance \/r^/{\ — XBj) is the associated transverse distance between the struck parton (probed 
during the hard interaction) and the center of momentum of the spectators. That's why it can be 
interpreted as a typical spatial extension of partons in the proton. 

What we have learned so far with the present experimental situation is already very rich: slow 
partons (at low XBj) are located at the periphery of the proton whereas fast partons (at large XBj) 
make up the core of the proton (in its center). This last property is only an indirect observation from 
fits of form factor measurements |29j (see below for a short discussion). 

We need to get more information. How large can be the spread in space of slow partons? Could 
it be larger that 1 fm? Also, what is the spread for the large x (constituent) partons? Where is the 
transition between the large x partons and the peripheric partons? We need more experimental results 
and then more experiments with different setups to address these questions from all possible angles. 

For example, if we would observe a gradual increase of the t dependence of the GPD H{x,0,t) 
(quarks or gluons) when varying xbj from large to small values, it would mean exactly that quarks at 
large xbj come from the more localized valence core of the proton, while the small x Bj region receives 
contribution from the periphery or, in other words, from the wider meson cloud. This is a very nice 
perspective for the future to expect direct measurements of \/r^ from many experiments in the world. 

10.4 An elegant application 

Let us come back briefly to form factors and their essential role in the interplay between x and t 
kinematic variables. A complete analysis is presented by Diehl et al. in Ref. [31]. Indeed, it is clear 
that indirect information on impact parameter distributions can be obtained by using the sum rules 
presented in Eq. (1181) . which provides a natural link between the GPDs dependences in x and t. We 
can exemplify the structure of the link on the Dirac form factor for proton and neutron 

Ff{t) = l\x[^H:{x,t) - ^H^{x,t)] 

Fi"(t) = j\x[^Ht{x,t) - \H:{x,t)] 

where we have neglected the contribution from the s quarks. Note that only valence type distribu- 
tions appear in these relations, since the electromagnetic current is only sensitive to the difference of 
quark and antiquark distributions. Then, from an ansatz for the functional dependence of if^(x, 0,t) 
and measurements of the Dirac form factor Fi{t) (and F2{t)), a fit of some GPDs parameters can be 
performed [3T] . 

Obviously, the sensitivity of such a fit is governed by the parameters building the interplay of x and 
t dependences (for valence distributions), which is the purpose of this approach. In Fig. HT] the default 
results for GPDs as tomography plots in impact parameter space is illustrated for fixed longitudinal 
momentum fraction x [31] . 

This confirms the results on \/r^ discussed above: low xbj partons are located at the periphery of 
the proton whereas valence like partons make up the core of the proton (in its center) . 
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Figure 47: Tomography plots of ^^(x, Ft) (left) and rx) (right) in the transverse rf-r^ plane. 
Note that the scale of intensity for longitudinal momentum fraction x = O.Q differs from the one for 
X = 0.3 and x = 0.05 BH. 
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Figure 48: Skewing factor R = XmA {'~f*p —>■ '~fp)t=o/Xrn A {'~f*p '~f*p)t=o extracted from DVCS and 
DIS cross sections [32] The GPD model is also displayed and gives a good agreement of the data (full 
line). The forward ansatz model used at all values of Q^, fails to reproduce the total skewing 
effects generated by the QCD evolution (dashed line). 



11 Quantifying skewing effects on DVCS at low xbj 
11.1 DVCS in the context of GPDs 

After this short overview of GPDs physics, we understand clearly why DVCS is the typical (and cleanest) 
process to extract GPDs, or at least to extract informations on GPDs. Then, we can come back on the 
DVCS cross section measurements and their interpretation in terms of GPDs. In order to quantify the 
magnitude of skewing effects, and thus the impact of GPDs on the DVCS process (7*p ■yp), we need 
to derive for example the following ratio from measured cross sections |l26j 

R = ImA{-f*p -fp)t=o/ImA{-f*p j*p)t=o- 

In this expression, ImA{^*V ~^ lP)t=o{Q^, is the imaginary part of the DVCS process and is directly 
proportional to the GPDs. Also, the diagonal term ImA {'y*p 'y*p)t=o is directly proportinal to the 
total cross section. The ratio R is then is equivalent to the ratio of the GPDs to the PDFs. That's why 
its measurement can provide directly the impact of GPDs, when compared to pure PDFs predictions. 



11.2 Experimental results 

In Ref. [32], we have shown how to extract this ratio from the DVCS and DIS cross sections. Results 
are presented in Fig. HHl The typical values of R are found around 2, whereas in a model without 
skewing R would be equal to unity. Therefore, the present measurements confirm the large effect of 
skewing. 

Values of R are also compared with a GPDs model based on a forward ansatz at low scale {Qq = 1.3 
GeV) [33]. Namely, the singlet GPD is parametrized as follows: Hs{x,^) = Qs{x), where Qs{x) is 
the singlet PDFs and x and ^ are the variables used in the previous part for the definition of GPDs 
(see Fig. HHD- It does not mean that the GPD is taken to be exactly the PDF. Indeed, at a; = ^, 
we get Hs{^,C,) = QsiO = Qsi.^Bj/'^) and not Qsi^Bj)- In other words, in this forward ansatz 
parameterization of the GPDs, we simply consider that at a low scale Qo, we can forget the profile 
function and take directly the parameterization of the GPD from a PDF like form. The same is done 
with non-singlet and gluon distributions. 
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If the GPDs are parametrized in such a way at initial scale, then we have two possibilities. Either, 
we evolve the GPDs using skewed QCD evolution equations, which naturally generates the skewing 
dependence (in ^) along the evolution, or we forget about the skewed evolution and we consider 
only the standard QCD evolution equations like for PDFs [32l [33] . This corresponds to the two curves 
presented in Fig. HHl The full line represents the complete GPDs model, with skewed evolution equations 
and the dashed curve, labeled forward ansatz (all Q^), represents the case where initial distributions 
are evolved with standard QCD evolution equations. Then, Fig. HH] demonstrates that we need the full 
GPDs model to describe our data on DVCS cross sections (converted in R values). If we forget about 
the skewing generated during the QCD evolution, we miss the data by about 30%. This is clearly a 
deep impact of the skewing effects present in the data [321 133] . 

Another influence of GPDs that we can check on data concerns the t dependence. We have already 
shown that in the kinematic domain of HI and ZEUS measurements, DVCS cross section {da /dt) can be 
factorized and approximated to a good accuracy by an exponential form e"'*'*' , which implies a factorized 
dependence in e"''/^'*' for GPDs. However, we can think of taking into account a non-factorized form 
in \x\~" as well. With the small a' value determined previously, we know that this term can only be 
small (negligible) correction to the dominant (factorized) t dependence in e"''/^'*' for GPDs. 

12 On the way of mapping out the GPDs 

12.1 Prospects for the COMPASS experiment at CERN 

As we have shown, the mapping of the GPDs is certainly a difficult work due to the flexibility of these 
functions. However, we have already illustrated some elements that can be constrained with the present 
DVCS data at low xsj- Concerning the t dependence of the GPDs in this kinematic domain, we have 
shown that the impact of a potential non-factorized term in is small, due to the small value of 

a' observed at low XBj- 

This is one important element of the experimental project to measure DVCS at COMPASS in the 
future, as we need to check this kind of effects at larger xsj- DVCS at COMPASS (located at CERN) 
can be measured with muon beams on fixed target, /xp — > /ip7. If the muon energy is large enough, 
for example = 190 GeV, DVCS dominates over the BH contribution (as for HI and ZEUS) so that 
DVCS cross section can be measured directly. 

At smaller lepton energy, = 100 GeV, the DVCS signal is not dominant and can not be measured 
directly. Then, we need to use the property that DVCS and BH, having identical final state, can interfere. 
When the DVCS cross section itself can not be measured, the interference can be observed. The strong 
interest is that the xbj kinematic domain of COMPASS follows the one of HI and ZEUS at larger xbj, 
with XBj ~ [0.05 — 0.15], thus much larger values than in the kinematic domain of HI and ZEUS. A 
project is ongoing to install a proton recoil detector in the COMPASS setup and then measure DVCS 
cross section or DVCS/BH interference p6]. Some tests have already been done to show the technical 
feasibility of the proposed experiment. 

Let us discuss how we can access an interference between DVCS and BH reactions. In fact, since 
these two processes have an identical final state, they can obviously interfere. The squared photon 
production amplitude is then given by 

1^1 ~ I^shI "I" I^dvcsI ~^ ^dvcs -^bh ~^ -^dvcs -^bhj (19) 

" V ' 

/ 

where Ag^ is the BH amplitude, A^iycs represents the DVCS amplitude and I denotes the interference 
term. 

For unpolarized proton target and lepton beam, the interference term can be written quite generally 
as a linear combination of harmonics of the azimuthal angle 0, which is the angle between the plane 
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Figure 49: Azimuthal distribution of the beam charge asymmetry measured at COMPASS at Ef^= 100 
GeV and \t\ < 0.6 GeV^ for 2 domains of xbj {xbj = 0.05 ± 0.02 and xbj = 0.10 ± 0.03) and 3 domains 
of Q2 (^Q2 ^ 2 ± 0.5 GeV^, = 4 ± o.5 GeV^ and = 6 ± 0.5 GeV^) obtained in 6 months of data 
taking with a global efficiency of 25% and with 2 ■ 10® ^ per SPS spill (P^+ = —0.8 and P^- = +0.8) 



containing the incoming and outgoing leptons and the plane defined by the virtual and real photons. 
In the leading twist approximation (at sufficiently high Q^), if only the first term in cos(j) and sin(^ are 
considered, it can be written as: 

/ oc —C [a cos (f) ReADvcs + bPi sin (p ImADvcs] (20) 

In this expression, C = ±1 is the lepton beam charge. Pi its longitudinal polarization and a and b are 
functions of the ratio of longitudinal to transverse virtual photon fiux [29] . 

At COMPASS, if we measure a beam charge asymmetry (BCA), the polarization of the muon beam 
fiips with the charge and so, the sin terms disappears. Then, the BCA reads 

da+/d(p + da /dcp \Advcs\ +\Abh\ 

Note that DVCS cross section measurements which are integrated over are not sensitive to the 
interference term (see Eq. |20|) . Simulations done for COMPASS [36] are shown in Fig. |l9]for BCA in a 
setup described in the legend of the figure. Two models of GPDs, with a factorized and non-factorized 
t dependence, are shown in Fig. HH and we can observe easily the great discrimination power offered by 
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Figure 50: Simulation of the azimuthal angular distribution of the beam charge asymmetry measurable 
at COMPASS at = 100 GeV. We present the projected values and error bars in the range \t\ < 0.6 
GeV^ for 2 values of xbj (0.05 and 0.1) at = 4 QeV^ (see Ref. |36]). The prediction of the GPD 
model with a non-factorized t dependence is shown (full line). The case of a factorized t dependence 
would lead to a prediction of the BCA compatible with zero and is not displayed. 



COMPASS, with the proton recoil detector fully operational [36]. Of course, the discrimination is large 
in Fig. HHldue to the fact that a' is taken to be large (a' ~ 0.8 GeV'^) in simulations. If it happens to 
be much smaller, as measured at low xbj by HI [26] (see previous section), both predictions for BCA 
in Fig. |l9] would be of similar shape, as both curves would converge to the factorized case. 

In Fig. \5U[ we compare predictions of the GPD model used in the previous section for HI data to 
simulations of the BCA extraction at COMPASS using a muon beam of 100 GeV [32]. We present the 
comparison for one value of (4 GeV^) and two values of xbj (0.05 and 0.1). When we compute the 
BCA in the factorized exponential t dependence approximation, we find values compatible with zero, 
which are not represented in Fig. [501 

Therefore, we display only the predictions of the model obtained in the non-factorized case using 
the same a' ~ 0.8 GeV~^ value than in Ref. [36]. Both the cos(0) and cos(2(/)) terms contribute to a 
significant level to the BCA at COMPASS, as illustrated in Fig. [501 We notice that our predictions do 
not match with the COMPASS simulation done with the model described in Ref. [SB]. This is another 
illustration of the large discriminative power of this observable on GPDs parameterizations, even for 
identical t dependence input. 



12.2 Recent results on azimuthal asymmetries at HERA 

At HERA, we have also samples of data with electron and positron beams. Therefore, it has been 
possible to extract the beam charge asymmetry, Ac = ^ former pioneering measurement 

of the BCA at HERMES [3l| is shown in Fig. [SH HERMES was a fixed target experiment located 
at DESY operating with the electrons or positrons beams of 27.6 GeV. Recent results from HERMES 
[3l] and HI [26] are presented in Fig. [52] and [531 which correspond to xbj ~ 0.1 for HERMES and 
XBj ~ 10~^ for HI. Note that for HI results, we have kept a different convention in the definition of 
than in fixed target experiments, namely (pni = vr — (Phermes- The advantage of the convention we 
have considered in HI is that, a positive pi (with BCA = picos0) means a positive real part of the 
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Figure 51: Beam charge asymmetry Ac from HERMES data as a function of \(f)\. Statistical uncertainties 
are shown. The sohd curve represents the four-parameter fit: (—0.011 ±0.019) + (0.060 ± 0.027) cos0 + 
(0.016 ± 0.026) cos 20 + (0.034 ± 0.027) cos 30. The dashed hne shows the pure cos0 dependence. 
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Figure 52: Moments of the Beam charge asymmetry from HERMES data. 
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Figure 53: Beam charge asymmetry as a function of from HI data ^2^. Note that the convention 

chosen in HI is different from HERMES (see text). The inner error bars represent the statistical errors, 
the outer error bars the statistical and systematic errors added in quadrature. A fit in picos0 is 
presented, where pi is a free parameter (full line), together with the GPDs model prediction (dashed 
line) . 

DVCS amplitude. 

Both experiments show that the present status of GPD models can correctly described the BCA 
measurements. In case of HI, factorized parameterizations of GPDs (in t) are the most simple choices 
compatible with measurements (see above), and for HERMES, the sensitivity of the hypothesis of the 
t-dependence is illustrated in Fig. [52l 

If we consider the overall description in Fig. [52], the factorized ansatz (without D-term) is favored by 
HERMES BCA measurements. The so-called D-term is part of some parameterizations of GPDs |34j . 
That's why BCA, which provides a sensitivity the real part of the DVCS amplitude, gets a sensitivity 
to this (unknown) term. In general also, the factorized ansatz is much more stable with respect to the 
D-term contribution, when compared to the non-factorized (Regge) ansatz. Indeed, the spread between 
Regge with/without D-term predictions is huge, whereas the D-term has only a small impact on the 
factorized predictions. As the D-term is almost completely unknown, it is interesting to make choice of 
parameterizations (if possible) that can reduce their sensitivity to it. 

In Ref. [3l], it is mentioned that the Regge (without D-term) is favored, based on the observation 
of the t dependence. However, it is not that clear when considering all data points. 

In any case, the experimental results presented in Fig. [52] and [53] are the first obtained on BCA and 
then important pieces to provide constraints on the real part of the amplitude in future developments 
of GPDs phenomenology. A compilation of HI and HERMES results is presented in Fig. [JU 

12.3 Experimental analysis of dispersion relations 

A specific analysis has been done in the HI experiment concerning the real part of the DVCS amplitude 
|26j . From Eq. f[2T]) and measurements of BCA and DVCS cross section, it is possible to extract the 
ratio of the real to imaginary parts of the DVCS amplitude 



This ratio is a key quantity which can also be derived through a dispersion relation, which takes a 
simple form in the high energy limit. Indeed, at low xbj, when the W dependence of the DVCS cross 



p = ReADvcsf^^^Dvcs- 
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Figure 54: The beam charge asymmetry for HI and HERMES is fitted by a the dominant harmonic 
in picos0, where pi is a free parameter (see text). Results for pi are presented for both experiments 
as a function of xbj- This plot can be interpreted as a reflection of the different values of a' in both 
kinematic ranges, negligible value in the HI domain and large value for HERMES. We understand 
immediately the great interest of COMPASS kinematical domain, lying at intermediate xbj- 



section is dominated by a single term in (with 6 > 0.3), the dispersion relation can be written as 



ReAovcs /'i-'^ADvcs = tan 



(22) 



where (5(Q^) is the power governing the W dependence of the DVCS cross section at a given Q^. As 
we have measured 6 independently from DVCS cross sections only (see previous section), we can 
compute this ratio, with the very reasonable assumption that the dispersion relation are correct. We 
obtain: p = 0.25 ± 0.06. To be compared to the value extracted from BCA measurement and the 
subsequent extraction of pi, which gives p = 0.23 ± 0.08. Both values are found in good agreement. 

After this brief discussion, we can also understand simply how the sensitivity of the beam charge 
asymmetry observable is built with a'. The BCA is proportional to the ratio of real to imaginary part 
of the DVCS amplitude and this ratio can be expressed with respect to t as 

p[t] = ReADvcsf^^ADvcsit] ^ ReADvcs/'i-TxiADvcs[0]il - a'|t|). 

Then, trivially, for small values of a' at low |t| values, we do not expect much sensitivity (on a') of this 
ratio and thus of the BCA. This is what is illustrated for HERMES results in Fig. 



12.4 Jefferson Laboratory experiments 

Regarding the kinematic coverage of fixed-target experiments (see Fig. [55]) . the Jefferson Lab (JLab) 
experiments play a major role in the field, exploring the large xbj and low kinematic domain. JLab 
experiments, colliding an electron beam in the energy range of 6 GeV on a fixed target, can measure 
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beam spin or target spin asymmetries [35] and then access directly the imaginary part of the DVCS 
amphtude in the valence domain. 

Of course, we can not exclude a priori that higher twists effects would completely spoil any pertur- 
bative treatment of the experimental results in this area. Below, we describe briefly few characteristic 
measurements at JLab related to GPDs physics. 

First, let us recall that in this kinematic domain, the BH cross section is completely dominating 
the ep ep7 cross section. Then, the DVCS signal can hardly be observed and only the BH/DVCS 
interference can be accessed through different observables with different sensitivities to GPDs [55] . 

An important recent result has been obtained by the Hall A E-00-110 experiment [3S], which demon- 
strates that measurements at JLab are dominated by leading twists contributions. It is shown in Fig. 
I56l From the observed scaling of the imaginary part of the DVCS amplitude (see Fig. [561) this result 
provides an indication that the measurement of the imaginary part of the DVCS amplitude follows a 
typical Bjorken scaling, observed over the range covered by the experiment. Which means that 
leading twists terms are likely to dominate. Of course, the range in accessible is not large but the 
the high precision of the data makes this last statement quite reasonable. An upgrade at larger energies 
of the lepton beam is obviously an important issue to get a sensitivity to higher values (and larger 
W). 




Figure 55: Kinematic coverage for fixed-target experiments: (i) Compass at 190 GeV; (ii) Hermes 
at 27.6 GeV, dotted line for existing data (< 2005), solid line for future (2005-2007) data with an 
integrated luminosity higher by about one order of magnitude; (iii) JLab experiments at 6 GeV, and 
at 11 GeV (after upgrade). 

Apart from DVCS/BH interference measurements, a separation between BH and DVCS processes 
has been obtained with the Hall A E-00-110 experiment. The measurement of the 4-fold (polarized and 
unpolarized) differential cross sections ^^^^^2^^^^ (for the real photon production process) has been done 
for three values of Q^{in the kinematic domain W ^ 2 GeV and x > 0.1). Results are shown in Fig. [57] 
for < >=2.3 GeV^ [35|. The particular shape in of the unpolarized cross section (Fig. [STj) is typical 
of the BH process. The dot-dot-dashed curve in Fig. [57] shows its precise shape and contribution. It 
can be seen that it dominates most of the cross sections and, only around $= 180°, there is a large 
discrepancy (a factor ?a 2) between the BH and the data which could be attributed to the DVCS process 
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Figure 56: dependence of imaginary part of the DVCS amplitude (left). We observe a scaling over 
the range in covered by the analysis. 




Figure 57: Difference of (beam) polarized cross sections for DVCS on the proton, as a function of 
the $ angle, measured by the JLab Hall A collaboration. The average kinematics is < xsj >=0.36, 
< >=2.3 GeV^ and < —t >=0.28 GeV^. The figure on the bottom shows the total (i.e unpolarized) 
cross section as a function of The BH contribution is represented by the dot-dot-dashed curve. 
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itself. It opens a possibility in a future analysis to extract directly a DVCS signal, which would be a 
first time measurement in this kinematic range. 




^ (deg) 

Figure 58: Left: kinematic coverage and binning in the {xb, Q"^) space. Right: Beam spin asymmetry 
A{(j)) for 2 of the 62 (xb, Q^, t) bins, corresponding to {xb) = 0.249, (Q^) = 1.95 GeV^, and two values 
of (t). The long-dashed curves correspond to fits with A = ^ f ^ • The dashed curves correspond to 
a Regge calculation. GPDs calculations are also shown as full lines. 



Let us present a final measurement from the (JLab) Hall B E-00-113 experiment, concerning beam 
spin asymmetries (BSAs) [35j, which shows (again) clearly the interest for an upgrade at larger energies. 
Results are presented in Fig. [5H] with GPDs or Regge (non-perturbative) models. The asymmetries are 
fitted according to the relation 

a sin 6 , , 

A = — (23) 

1 + c cos (p + a cos 2(p 

and extracted values of a are displayed in Fig. [591 As can be seen in Fig. [591 the discrimination of Regge 
(soft) or GPDs (hard) approaches is not conclusive from the present data. Therefore, the upgrade with 
12 GeV electrons is very interesting to address this separation between soft and hard physics at JLab. 



12.5 Experimental prospects on the orbital angular momentum of partons 

A final comment is in order concerning the measurement of asymmetries (from DVCS/BH interference) 
in fixed target experiments. Experiments at JLab and data collected by HERMES allow to determine 
transverse target-spin asymmetries, by controlling the polarization of the target. This would be also 
a possibility of the future COMPASS project described above. Experimentally, we need to introduce 
another azimuthal angle (ps to characterize completely the events measured in such configurations, 
where (ps represents the direction of the spin of the target with respect to the plane of the leptons 
(incident and scattered). 

The great interest is then that the cos(j) moment of the asymmetry da{(j),(j)s) — dcr{(f),(j)s + tt) is 
proportional to the imaginary part of GPDs of types H and E. Remind the short note we have written 
in the last section: the contribution of GPDs of type E are damped by kinematic factors of orders 
\t\/Mp in all the observables we have discussed till now. This is not the case for transverse target-spin 
asymmetries. 

Therefore, these measurements are particularly interesting in the quest for GPDs. The strong interest 
in determining GPDs of type E is that these functions appear in a fundamental relation between GPDs 
and angular momenta of partons. Indeed, GPDs have been shown to be related directly to the total 
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Figure 59: Extracted values a of the DVCS-BH interference as measured in a spin asymmetry, with 
A = T: — "l™!^ — FTT. They are shown as a function of —t for different bins in (xr, Q^). The sohd curves 

l+ccos (p+rfcos 2(p -J \ Dl / 

represent a GPD model and the dashed curves a Regge approach. 



angular momenta carried by partons in the nucleon, via the Ji relation [29] 

1 

-j^dXX (HgiX, e, t) + Eg{X, e, t)) = Jg. (24) 

As GPDs of type E are essentially unknown apart from basic sum rules, any improvement of their 
knowledge is essential. From Eq. fl2^ . it is clear that we could access directly to the orbital momentum 
of quarks if we had a good knowledge of GPDs H and E. Indeed, Jq is the sum of the longitudinal 
angular momenta of quarks and their orbital angular momenta. The first one is relatively well known 
through global fits of polarized structure functions. It follows that a determination of Jq can provide an 
estimate of the orbital part of its expression. In Ji relation (Eq. ( 12^ ). the function H is not a problem 
as we can take its limit at ^ = 0, where H merges with the PDFs, which are well known. But we need 
definitely to get a better understanding of E. 

First measurements of transverse target-spin asymmetries have been realized at JLab [35| and HER- 
MES We present results obtained by HERMES ^341 in Fig. [HDl The typical sensitivity to hypothesis 
on Jq values is also illustrated in Fig. EHl with the reserve that in this analysis, the observed sensitivity 
to Jq is model dependent. It is already a first step, very challenging from the experimental side. Cer- 
tainly, global fits of GPDs (if possible) would give a much more solid (less model dependent) sensitivity 
to Jq (see next section). 

12.6 A few comments on the Ji relation 

In order to give more intuitive content to the Ji relation (^^, we can comment further its dependence 
in the function E. From our short presentation of GPDs, we know that functions of type E are related 
to matrix elements of the form {p', s'\0\p, s) for s ^ s', which means helicity flip at the proton vertex 
(s 7^ s'). That's why their contribution vanish in standard DIS or in processes where t tends to zero. 
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Figure 60: Target-spin asymmetry amplitudes describing the dependence of the squared DVCS amph- 
tude (circles, Aut,dvcs) and the interference term (squares, Aut,i) on the transverse target polarisation. 
In the notations, U refers to Unpolarized beam and T to Transversely polarized target. The circles 
(squares) are shifted right (left) for visibility. The curves are predictions of a GPD model with three 
different values for the w-quark total angular momentum and fixed d-quark total angular momentum 
Jrf = (see [M])- This is a first important (model dependent) check of the sensitivity these data to the 
Ji relation. 



More generally, their contribution would vanish if the proton had only configurations where helicities 
of the partons add up to the helicity of the proton. In practice, this is not the case due to angular 
momentum of partons. This is what is reflected in a very condensed way in the Ji relation (Eq. (!24|) ). 

Then, we get the intuitive interpretation of this formula: it connects E with the angular momentum 
of quarks in the proton. A similar relation holds for gluons [29], linking Jg to Hg and Eg and both 
formulae, for quarks and gluons, add up to build the proton spin 

Jg + Jg = 1/2. 

This last equality must be put in perspective with the asymptotic limits for Jq and Jg at large scale 
Q^, which read Jq \ le+snj ^^'^ "^s ~^ \ i6+3nf ' where % is the number of active flavors of quarks at 
that scale (typically = 5 at large scale Q"^) [29] . 

In words, half of the angular momentum of the proton is carried by gluons (asymptotically). It is not 
trivial to make quantitative estimates at medium scales, but it is a clear indication that orbital angular 
momentum plays a major role in building the angular momentum of the proton. It implies that all 
experimental physics issues that intend to access directly or indirectly to GPDs of type E are essential 
in the understanding of the proton structure, beyond what is relatively well known concerning its 
longitudinal momentum structure in XBj- And that's also why first transverse target-spin asymmetries 
(which can provide the best sensitivity to E) are so important and the fact that such measurements 
have already been done is promising for the future. 

Clearly, we understand at this level the major interest of GPDs and we get a better intuition on 
their physics content. They simultaneously probe the transverse and the longitudinal distribution of 
quarks and gluons in a hadron state and the possibility to flip helicity in GPDs makes these functions 
sensitive to orbital angular momentum in an essential way. This is possible because they generalize the 
purely collinear kinematics describing the familiar twist-two quantities of the parton model. This is 
obviously illustrating a fundamental feature of non-forward exclusive processes. 
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Figure 61: (a) Global GPD fits at t = —0.3 GeV^ (thick) and t = (thin) are displayed as dashed and 
sohd lines, respectively, (b) Prediction of the BCA for COMPASS kinematics (S^, = 160 GeV, = 2 GeV^, 
t = -0.2 GeV2) [301. 



12.7 Towards global fits in the GPDs context 

A direct continuation of the analysis exposed in the previous section is to perform a global fit of all 
previous experimental results. In the same spirit as it is done for global QCD fits of proton structure 
function data, obtained in DIS scattering, a global fit can be done of observables measured for exclusive 
processes, like exclusive real gamma production. Instead of defining initial conditions on PDFs (at a 
low scale Qq), initial conditions on GPDs must be assumed. An important step in this direction is 
presented in Ref . [30] . 

A typical result derived in this work is displayed in Fig. [6ll The GPD of type H is shown for two 
values of t (see FigEU-left-) and the influence of JLab results is illustrated on the prediction of the BCA 
in the COMPASS kinematics (see Fig|6T]-right-). 

From these global GPDs fits [l30J, the impact parameter space distribution can be extracted with 

p(rT, X, Q2) = 1^ ^ "-Ll \ ^ (25) 

Results obtained in Ref. [30] confirms what we have already discussed in previous sections [3Tj. This 
theoretical framework to analyze GPDs is a promising trend for the future, in parallel to the production 
of new experimental measurements. 



13 Outlook 

We have reviewed the most recent experimental results from hard diffractive scattering at HERA and 
Tevatron. We have shown that many aspects of diffraction in ep collisions can be successfully described in 
QCD if a hard scale is present. A key to this success are factorization theorems, which render parts of the 
dynamics accessible to calculation in perturbation theory. The remaining non-perturbative quantities, 
namely diffractive PDFs and generalized parton distributions, can be extracted from measurements 
and contain specific information about small-x^j partons in the proton that can only be obtained in 
diffractive processes. To describe hard diffractive hadron-hadron collisions is more challenging since 
factorization is broken by re-scattering between spectator partons. These re-scattering effects are of 
interest in their own right because of their intimate relation with multiple scattering effects, which at 
LHC energies are expected to be crucial for understanding the structure of events in hard collisions. 
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A combination of data on inclusive and diffractive ep scattering hints at the onset of parton saturation 
at HERA, and the phenomenology developed there is a helpful step towards understanding high-density 
effects in hadron-hadron collisions. In this respect, we have discussed a very important aspect that 
makes diffraction in DIS so interesting at low xbj- Its interpretation in the dipole formalism and its 
connection to saturation effects. Indeed, diffraction in DIS has appeared as a well suited process to 
analyze saturation effects at large gluon density in the proton. In the dipole model, it takes a simple and 
luminous form, with the introduction of the so-called saturation scale Qs- Diffraction is then dominated 
by dipoles of size r ~ 1/Qs- In particular, it provides a simple explanation of the Constance of the ratio 
of diffractive to total cross sections as a function of W (at fixed values). 

Then, exclusive processes in DIS, like VMs production or DVCS, have appeared as key reactions 
to trigger the generic mechanism of diffractive scattering. Decisive measurements have been performed 
recently, in particular concerning dependences of exclusive processes cross section within the momen- 
tum exchange (squared) at the proton vertex, t. This allows to extract first experimental features 
concerning proton tomography, on how partons are localized in the proton. It provides a completely 
new information on the spatial extension of partons inside the proton (or more generally hadrons), as 
well as on the correlations of longitudinal momenta. A unified picture of this physics is encoded in the 
GPDs formalism. We have shown that Jefferson laboratory experiments or prospects at COMPASS 
are essential, to gain relevant information on GPDs. Of course, we do not forget that the dependence 
of GPDs on three kinematical variables, and the number of distributions describing different helicity 
combinations present a considerable complexity. In a sense this is the price to pay for the amount of 
physics information encoded in these quantities. It is however crucial to realize that for many impor- 
tant aspects we need not fully disentangle this complexity. The relation of longitudinal and transverse 
structure of partons in a nucleon, or of nucleons in a nucleus, can be studied quantitatively from the 
distribution in the two external kinematical variables Xsj and t. 
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Figure 62: Quark spin and orbital angular momentum contributions to the spin of the nucleon for up 
and down quarks. Squares and triangles denote AS" and AS*^ respectively, and diamonds and circles 
denote L" and L'^ respectively [37] . 



To conclude, we can illustrate these issues with results from lattice QCD [37J. In Ref. [37], lattice 
QCD calculations are performed. They show two remarkable features of the quark contributions to the 
nucleon spin. The first is that the magnitude of the orbital angular momentum contributions of the up 
and down quarks, and L'^, are separately quite substantial, and yet they cancel nearly completely 
(see Fig. [62l) . The second is the close cancellation between the orbital and spin contributions of the d 
quarks, L'^ and AE'^/2. Of course, we can not take these results as granted but calculations are solid. 
It would be obviously valuable to understand the physical origin of both features, with more data. 
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